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ABSTRACT 

WMAP data when combined with ancillary data on free-free, synchrotron and dust 
allow an improved understanding of the spectrum of emission from each of these 
components. Here we examine the sky variation at intermediate and high latitudes 
using a cross-correlation technique. In particular, we compare the observed emission 
in several large partitions of the sky plus 33 selected sky regions to three "standard" 
templates. The regions are selected using a criterion based on the morphology of these 
template maps. 

The synchrotron emission shows evidence of steepening between GHz frequencies 
and the WMAP bands. There are indications of spectral index variations across the 
sky but the current data are not precise enough to accurately quantify this from 
region-to- region. 

The emission correlated with the Hq, template shows clear evidence of deviation 
from a free-free spectrum. The emission can be decomposed into a contribution from 
both free-free and spinning dust in the warm ionised medium of the Galaxy. The 
derived free-free emissivity corresponds to a mean electron temperature of ~ 6000 K, 
although the value depends critically on the impact of dust absorption on the Hq 
intensity. The WIM spinning dust emission has a peak emission in intensity in the 
range 40-50 GHz. 

The anomalous microwave emission associated with dust is detected at high signif- 
icance in most of the 33 fields studied. The anomalous emission correlates well with the 
Finkbeiner et al. (1999) model 8 predictions (FDS8) at 94 GHz, and is well described 
globally by a power-law emission model with an effective spectral index between 20 
and 60 GHz of /3 w —2.7. It is clear that attempts to explain the emission by spinning 
dust models require multiple components, which presumably relates to a complex mix 
of emission regions along a given line-of-sight. An enhancement of the thermal dust 
contribution over the FDS8 predictions by a factor 1.2 is required with such models. 
Furthermore, the emissivity varies by a factor of ^ 50% from cloud to cloud relative 
to the mean. 

The significance of these results for the correction of CMB data for Galactic fore- 
ground emission is discussed. 

Key words: cosmology: observations - cosmic microwave background - radio contin- 
uum: ISM ~ diffuse radiation - radiation mechanisms: general 
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1 INTRODUCTION 

A major goal of observational cosmology is to determine those parameters that describe the Universe. Observations of the 
Cosmic Microwave Background (CMB) at frequencies in the range 20 - 200 GHz provide unique data to achieve this by 
establishing the statistical properties of temperature (and polarisation) measurements. However, an impediment to such 
studies arises due to foreground emission in our own Galaxy from at least three sources - synchrotron, free-free and thermal 
dust emission. As CMB studies move to ever higher precision it is essential to determine the properties of these components 
to similarly high accuracy. Indeed, although the combined foreground level reaches a minimum in this frequency range 
{f ~ 70 GHz), it remains the dominant signal over large fractions of the sky. Of particular relevance to this discussion is the 
fact that each of the foreground components has a spectral index that varies from one line of sight to another so using a single 
spectral index can lead to significant uncertainties in the corrections required. It is therefore essential to study the nature of 

the Galactic signal at microwave wavelengths in its own right. 

All-sky observations by the Wilkinson Microwave Anisotropy Probe ( WMAP, iBennett et a l."2003a'') at the 5 frequencies 
of 23, 33, 41, 61 and 94 GHz can provide the basis for improving our understanding of local foregrounds. By comparing these 
maps with templates for synchrotron, free-free and dust emission, made at frequencies where the specific emission mechanisms 
dominate, it is possible to clarify important properties of the emission. Indeed, new insights into the nature of the Galactic 
diffuse emission have arisen from studies o f the WMAP data , including the detection of several unexpected new c ontributions. 

Anomalous dust-correlated emission (|Leitch et al.lll99 7i) was originally observed in the CO-B-E-DMR data (Kogut et cdj 

Il99d ) but was thought to be due to free- free emission. Draine & Lazarian (1998a,b) shifted attention to the dust itself as the 
source of emission through a mechanism now referred to as "spinning dust" , or dipole emission from very rapidly spinning 
grains. A rea nalysis of the in termediate and high Galactic latitude data taken by COS-E-DMR and supplemented by 19 GHz 
observations jBandav et al.l 2003) led to evidence for dust at intermediate Galactic latitudes emitting a spectrum consistent 
with the the form expected for spinning dust, specifically a hint of a turnover at frequencies below ~ 20 GHz. However, 
study of the WMAP data has allowed f urther refinement of our unde rstanding of the emission. Cross-correlation of the K- 
band data with observations at 15 GHz (|de Oliveira-Costa et al.l 2004) ag ain indicated a plateau or downturn in foreground 



emission inconsistent with a free- free or synchrotron origin. 'Lagach^ ( 2003 ) compared the WMAP data to HI column density 



measurements, revealing an increase in emission with decreasing density, suggesting that the anomalous emiss ion is connected 
to small, transient heated grains. In addition, a series of papers (|Finkbeineill20(Ml : lOobler fc Finkbeineilliooial lbl) have strongly 
confirmed the presence of anomalous dust emission in the WMAP data, and claim to have found evidence of such a component 
from the diffuse warm ionised medium (WIM) of the Galaxy. Specifically, the correlation with a H^ template, commonly utilised 
as a proxy for the free-free emission, is not consistent with the spectrum expected for ionised gas, and a broad bump is seen 
peaking towards ~ 40 GHz. Subsequently, Doblcr, Drainc & Finkbcincr (2009) have attempted to constrain specific spinning 
dust parameters such as the density and typical electric dipole moment of the grains. Recently, iPeeJ l|201lf) have shown that 
the K-band dust-correlated component is not strongly affected by the inclusion of a 2.3 GHz synchrotron template, reducing 
the possibility of a flat-spectrum synchrotron component. 

Finally, the so-cal led VFMAP-ha z e was identified bv lFinkbeiiieil (|200i) although it was alrea dy clear l y apparent in the 
foreground residuals of lBennett et al.l (|2003bh and subsequently in the SMICA analysis of .Patanchon et al.l (|2005l ) . The initial 
physical interpretation of the haze was that it was associated with free-free emission from hot gas in excess of lO'^ K, but this 
was refuted on the basis of the lack of associated X-ray emission. The detection of the haze relies upon the use of standard 
templates to remove known foreground emission utilising spatially independent spectra over the entire (high-latitude) sky. It 
has been argued that such a crude approxim ation to the true b ehaviour of the foreground emission at microwave frequencies 
may well lead to unphysical results. Indeed, iGold et al.l l|201ll l find that a spatial variation of spectral index of order 0.25 
between 408 MHz and K-band is sufficient to reproduce the haze amplitude. Furthermore, a corresponding polarised signal 
was found to be absent from the WMAP data. However, the lack of polarised emission can be explained by the entanglement 
of the Galactic magnetic field towards the Galactic centre, leading to a lower polarisation fraction there as compared to the 
outer Galaxy. Nevertheless, the microwave haze remains an active area of resea rch, particularly given its possible association 
with a gamma-ray counterpart observed in the Fermi data (|Dobler et al.|[201oh . 

In this paper, we characterise the spatial variation in the foreground emission in terms of the various emissio n mechanisms 
noted above by introducing a new partition of the sky into morphologically selected regions. In previous work jPavies et al.l 
I2O.O6, hereafter D06), our approach was to identify regions away from the Galactic plane which were expected to be dominant 
in one of the three foreground components, free-free, synchrotron or dust and to derive the spectrum for each component. 
Five regions covering angular scales of 3deg to 20deg were chosen for each component, based on foreground template maps, 
making 15 in all. Here, we generalise this approach and introduce an algorithm to define a set of 35 regions for further study. 
The selection is intended to minimise the potential cross-talk between the various physical components and to select regions 
over which the spectral behaviour is uniform thus supporting the use of a template-based comparison. We use a well-known 
and understood cross-correlation technique to undertake the analysis. 

The paper is organised as follows. Section [5] describes the WMAP data and foreground templates used in this analysis 
while Section [2. 31 defines the regions of interest for investigation. The methodology of the cross-correlation analysis is outlined 
in Section [3] and the corresponding results are presented in Section [l] Model-dependent spectral fits for each component 
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are considered in Section [5] and overall conclusions given in Section |6] Appendix |X] discusses in detail issues related to the 
Ha template used in the analysis, Appendix [B] defines the detailed method for partitioning the sky, and finally Appendix [Cl 
tabulates all of the template-fit coefficients for all the regions analysed. 



2 DATA USED IN THE ANALYSIS 
2.1 WMAP data 

We use the WMAP seven-year data (jjarosik et al.ll201ll ) provided in the HEALPi>Q pixelisation scheme l|G6rski e t al."2005l) 
with associated resolution Nside=^^'2 that can be obtained from the Legacy Archive for Microwave Background Data Analysis 
(LAMBDA) websitel!]. The WMAP satellite has 10 so-called differencing assemblies (DAs) distributed over five frequencies 
from ~23 GHz (K-band) to ^ 94 GHz (W-band), and with frequency-dependent resolutions ranging from approximately 0°.93 
up to 0°.23. The K- and Ka-bands have only one DA each, Q and V band have two each, while W-band has four. 

Multiple DAs at each frequency for the Q-, V- and W-bands are combined using simple averaging to generate a single map 
per frequency band. The data are then smoothed to a common resolution of 3° (after deconvolving the effective azimuthally 
symmetric beam response for each map) and degraded to HEALPix Nsid<i=6^. We perform our analysis at this resolution, 
rather than the more typical 1° studies, to account for the full covariance properties of the signal (see section We have 
also identified a problem with using the H^ data at higher resolution, as discussed further in section 12.21 and Appendix [X] 
Finally, we convert the data to brightness (antenna) temperature units from thermodynamic temperature since this is more 
appropriate for studying the spectral dependences of foregrounds. 



2.2 Foreground templates 

Each of the templ ates u sed here has been discussed at length in many prior studies. Thus we provide only a brief review and 
refer the reader to Ido^ and references therein. 

Synchrotron Template: The synchrotron emission arises mainly due to the acceleration of relativistic cosmic ray electrons in 
the Galactic magnetic field. Thus, the brightness temperature of the synchrotron spectrum depends on the energy spectrum 
of the cosmic ray electrons and the strength of the magnetic field. For an electron population with energy distribution N{E) 
given by a power law, N{E)dE oc E~^dE , the brightness temperature of the ensemble synchrotron spectrum is also a power 
law given by, r(i/) cx v^^' , where fis is related to spectral index of the energy spectrum 3 by the relation Ps — {5 + 3)/2. At 
very low frequency (< 1 GHz), the observed sky signal is dominated by the synchrotron emission from our Galaxy, and is least 
contaminated by free-free emission, at least away from the Galactic plane. The spectral behaviour changes with frequency as 
a result of the details of cosmic-ray electr on propagation. Ind eed, the mean spectral index changes typically from jSs — —2.55 
to —2.8 at 38 and 800 MHz respectively fLawson et al."l987'). Moreover, 13s varies across the sky and a range of values from 
-2.3 to -3.0 has been deterni ined by Reich & Reich (1988) between 408 and 1420 MHz data. The 408 MHz radio continuum 
survey of iHaslam et al.l l|l982t) provides a full sky template for synchrotron studies at an angular resolution close to 1°. For 
our purpose, we use the version of the data provided by the LAMBDA site. 

Free-Free Template: Free-free emission arises in regions of ionised hydrogen and is produced by free electrons scattering from 
ions without capture. The intensity spectrum of the free-free emission depends on the electron temperature and emission 
measure (EM) which is related to the number density of electrons (rie) along a given line of sight as EM= J n^dl. The optical 
Ha recombination line results from the capture of free electrons by a proton imcleus and is therefore also related to the EM. 
A high resolution full sky map of H^ emission can then be used as a good tracer of the free-free continuum emission at 
radio wavelengths. There is a well-defined relationship between the H^ intensity and radio- continuum bri ghtness temperature 
with a strong dependence on frequency and modest dependence on electron temperature Tp . iDickinson. Da vics & Davis (200^ 



herea fter DDD) generated a ful l-sky Hg map as a co mposite of WHAM data in the northern sky (jRevnolds et a l. 1998; HaffnerJ 
I1999D and the SHASSA survey l|Gaustad et al.ll200ll ) in the south. We refer the reader to their paper for more complete details. 

A significant uncertainty when using the H^ te mplate is the absorption of Hg b y fore ground dust. The absorption can be 
estimated by using the dust column density maps of ISchlegel. Finkbeiner fc David (Il998l . hereafter SFD), and the parameter 
fd corresponding to the fraction of dust in front of the Hg that causes the absorption. |DDD showed that for Galactic longitudes 
in the range I — 30° - 60° and latitudes | & | = 5° - 15°, fd ~ 0.3. However, for local high latitude regions such as Orion 
and the Gum nebula there is little or no absorption by dust. This was also the interpretation favoured for high latitude He 
template fits to the COBE-T>MR data in lBandav et al.l l|2003f ). We adopt a default value fd = in this paper, but consider 
the impact of varying the absorption parameter on our results. We note that there is an effective degeneracy between fd 
and electron temperature Te for the interpretation of the fits between the WMAP data and the He template, and that the 
majority of such analyses assume that the dust is coext ensive with the H g emission, i.e. fd ~ 0.5. 

An alternative He template has been assembled bv lFinkbeineJ (|2003l ) using additional data that provides higher angular 



^ http://healpix.jpl.nasa.gov 

^ http://lambda.gsfc.nasa.gov 
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resolution on limited areas of the sky. We have examined and compared these two templates in detail and discuss the 
comparison in Appendix |X] In summary, we find that there are inconsistencies in the template fit results obtained with the 
two templates that are resolved only by smoothing the data to a resolution lower than 1°. For our main results we choose a 
resolution of 3°. Such differences were already visible in the results of iDOS . These may be connected with a subtle interplay 
between artefacts in the templates and the template fitting method we utili se that leads to unstable estimates of the 
amplitude of the correlated emission at a given frequency. We use exclusively the IDDDI template fits in this work since the 
template was explicitly constructed to have uniform resolution. 

Dust Template: Thermal dust emission is the dominant foreground at frequencies of 100 - 1000 GHz. Its emissivity is generally 
modelled by a modified blackbody spectrum, I (i/) oc v"''^ B^(Td), where a d is th e emissivity index and B,j{Td) is the 
blackbody emissivity at a dust temperature Td- iFinkb einer. Davis fc Schlegell l| 19991 . hereafter FDS) predicted the thermal 
dust contribution at microwave frequencies from a series of models based on the COBi?-DIRBE 100 and 240 /im maps tied to 
COBS-FIRAS spectral data. The preferred model 8 (FDS8) has an effective power-law spectral index in antenna temperature 
of Pd ^ +1.55 over the WMAP frequencies. We use the FDS8 predicted emission at 94 GHz as our reference template for dust 
emission, and demonstrate that it also traces the anomalous component that predominates in the 10-100 GHz frequency range. 
Fits of the template to the WMAP data will help to constrain the spectral dependence of the dust-correlated foregrounds, 
and allow comparison to spinning dust models . Neve rtheless, we note that an alternative for the spinning dust component 
may ultimately be require d. Ind eed, iFinkbeinerl (|2004l ) proposed that an a nomalous dust tem plate could be better constructed 
using FDSSxTd^, however iDOfil find little evidence for this. More recently. lYsard. Miville-De sc henes fc Verstraete (^2010) have 
found improved correlation between the AME and the 12 /im brightness divided by the intensity of the interstellar radiation 
field in 27 fields of area 5° squared. However, such a template is not available at present for the full-sky given the difficulty 
in adequately cleaning the 12 fj,m. data from residual zodiacal emission. 

2.3 Masks for Global fits 

In order to examine the largest scale properties of the foregrounds, we have performed fits over large sky areas defined by 
two basic masks (EBVlfl and KQ85) and the intersection of the former with the Northern and Southern hemispheres defined 
in various reference frames (Galactic and ecliptic). We have also tested the sensitivity of the Northern hemisphere fits to 
the inclusion or otherwise of the North Polar Spur (NFS) region defined with reference to the 408 MHz radio survey. A 
representative subset of the effective masks are shown in Fig. [T] We collectively refer to analyses made using these masks as 
global fits. 

The KQ85 mask was developed by the WMAP team for the primary purpose of cosmological analysis. Since we are 
interested in investigating the properties of the foreground sky, our default EBV base mask is not as conservative as the KQ85 
mask that may exclude some areas of interest in our analysis. 

We define the EBV mask by rejecting that fraction of the sky for which the dust absorption in the H^ survey is unreliable. 
In practice, this correspo nds to the sky area where the absorption, as related to the dust optical depth maps from ISFd and 
the conversion factors in IDDDI exceeds 1 magnitude. This analysis is carried out at 3 degree FWHM resolution, then the 
extinction map is downgraded to a HEALPix pixel resolution Nside = 64. This is then merged with the downgraded WMAP7 
processing and point source masks in which any partially filled low resolution pixels are explicitly masked. Some additional 
pixels are also excluded corresponding to regions around the LMG, SMC and high-latitude clusters where some signal appears 
to leak outside the mask when applied to the WMAP7 K-band data smoothed to 3 degrees FWHM. 

For reference, we have included some fits based on the KQ85 mask. The Galactic plane part of the mask was smoothed 
to 3 degree resolution, then downgraded to Nside = 64 and those pixels with a value less than 0.95 were then rejected (set to 
0) and those above were accepted (set to 1). This was subsequently merged with the point source mask as used to create the 
EBV mask. 

Finally, the NFS may impact the synchrotron spectral indices determined in the Northern hemisphere. To test this, we 
apply a crude masking of the NFS. This is based on the observation that the 408 MHz map smoothed to 3 degrees resolution 
and divided by the cosecant of the Galactic latitude (an approximation to remove the path-length dependence of the emission 
within the Galaxy), shows an obvious enhancement towards the NFS compared to the rest of the sky. Thresholds are then 
applied to the map to eliminate data until there are no pixels in the North with values exceeding the maximum in the South 
(for Galactic latitudes less than —5 degrees). The resulting mask is certainly visually compelling as seen in Fig. [1] 

2.4 Regions Definition 

In this paper, we define a new set of regions generalising the approach from iDavies et all l|2006l ). In this previous analysis, 
regions were selected on the basis that one of the three foregrounds (free- free, dust or synchrotron emission as traced by 
standard templates) was dominant in each region. We extend this morphological definition in Appendix|Bj the result of which 
is the set of 35 candidate regions seen in Fig. [2l For analysis purposes, the regional mask is combined with the EBV mask, and 

^ The dust absorption correction is calibrated against an optical reddening law, or visual extinction, based on the (B - V) colours of a 
sample of galaxies, hence the name EBV for the mask. 
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Figure 1. A representative set of global masks used on this analysis. Masked areas are in blue. From left to right and from top to bottom 
: (1) EBV mask with point source mask included, (2) Galactic North Mask, (3) Ecliptic North Mask, (4) Galactic North Mask with 
North Polar Spur(NPS) removed (5) KQ85 mask with point source mask included, (6) Galactic South Mask (7) Ecliptic South Mask, 
(8) For reference, this shows the effective mask that would result if all of the 33 regions were combined and the EBV mask subsequently 
imposed. The individual regions are shown in Fig. [2] Note that this mask has not been used for analysis. 



this renders two candidate regions sufficiently small that we elect to omit them, leaving 33 regions for further study. These 
two regions can be seen in the figure but are not enumerated. Note that region 3 contains the NPS, but we do not attempt 
to suppress its impact by application of the NPS mask. In addition, our analysis is largely insensitive to the WMAP haze 
given that the sky areas where this is most prominent are not selected by our region definition scheme. This should be of 
little surprise given that the definition is based on the 408 MHz data in part, and haze emission is not seen therein. Various 
properties of the regions are specified in Table IBll 

We stress that we do not claim any absolute optimality of these regions for analysis. Indeed, they are by no means 
unique, and many alternate methods for regional definition are plausible. In fact, we recognise that, although we have applied 
a cosecant flattening to the data to allow their definition, such a flattening has not been applied in analysis and the flts may 
still be sensitive to the largest-scale Galactic emission as a function of latitude. Given the morphology of our regions, and the 
large range of latitudes that some cover (particularly towards the North Galactic pole) this could be a potential consequence of 
our approach. However, if there were a strong effect, it would be seen by comparing intermediate latitudes with high latitudes. 



3 METHODS 

The aim of this work is to cross-correlate datasets at different wavelengths in regions of the sky defined by the morphology 
of the templates described in § [2l The cross-correlation (C-C) method used here is a least-squares fit of one map to one or 
more templates. We perform the analysis at a resolution of 3° FWHM on sky maps d egrad ed to a HEALPix pixel resolution 
specified by Naide = 64 rather than the more typical 1° analysis seen, for example, in lDOGl This allows us to take account of 
the full covariance properties of the signal, and also provides a solution to some issues related to the the Hq, data at higher 
resolution, as discussed further in Appendix |X] In this way, we can include information about the CMB through its signal 
covariance rather than having to correct for it. 

It should be noted that in the high signal-to-noise regime probed by WMAP7 data, th e CMB covariance dominates the 
error budget in the fits. A different approach, as advocated for example bv lFinkbeinerl (l2004l ): lDobler fc Finkbeinerl (l2008al lbl). 
removes an estimate of the CMB signal from the data before cross-correlation. Such an internal linear combination (ILC) 
map typically corresponds to a linear combination of the data at the 5 WMAP frequencies. However, in principle this also 
correlates the fits between frequencies due to the common noise properties of the ILC map. More seriously, since the ILC 
contains foreground residuals, the subtraction changes the relative levels of foreground emissions at each frequency depending 
on the spectral characteristics of a given component. An app roach to rectify this has been specified in lDobler fc Finkbeinerl 
However, we prefer to retain the approach from iDOq in this analysis. 



3.1 Cross - Correlation Analysis 

The cross-correlation measure, a, between a data vector, d and a template vector t can be measured by minimising: 

X^ = {d- atf ■ MSl, ■ (d - at) = ■ M^^ • d (1) 

where Mgjv is the covariance matrix including both signal and noise for the template-corrected data vector d = d — at. 
Solving for a then becomes: 
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Figure 2. The set of regions defined in Appendix [B] 35 candidate regions have been identified, but two (here in dark blue and not 
enumerated) were considered to be too small after the EBV mask has also been applied for analysis purposes. 33 regions are analysed in 
this paper. 



(2) 



To compare multiple template components tj, e.g., different foregrounds, to a given dataset, the problem becomes a 
matrix equation. In the case where we have A'^ different foreground components, we end up with the simple system of linear 
equations Ax = b,where 



^kj — tfe 
Xk = 



MsN ■ d, 



(3) 



When only one template is present, this reduces to equation [5] above. 

The signal covariance is that for theoretical CMB anisotropies, Mfj = ^ X]fco(2^ ~^ ^)CiBj Bj^j^^Piifii ■ fij), where Bi 
is the beam window function for a Gaussian beam of 3° FWHM, and Bt^pi^ i s the HEALPix win dow function. The power 
spectrum, Ce, is taken from the WMAP best fit ACDM power law spectrum (jjarosik et al.|[201lh . The noise covariance is 
determined from the uncorrelated pixel noise as specified for each pixel in the WMAP data, and subsequently convolved 
and degraded appropriately to match the processing of the signal. The properties of the covariance matrix were also verified 
against noise simulations processed in the same way as the temperature maps. 

In this paper, we utilise three templates to describe the Galactic emission mechanisms, but also add a further template 
to describe monopole terms. Significant residual monopoles and dipoles may be present in the WMAP data, and can impact 
foreg round studies, particularly when based on parametric approaches (see for example lEriksen et al. I [2008; Dickinson et al.l 
|2009| ). Here, template fitting relies on the morphological content of the data, and the effective mean emission within a given 
region provides no relevant information. The inclusion of a single mon opole term can accoun t for offset contributions in all 
templates and the WMA P data in a way that does not bias the results (|Macellari et al.ll201ll ). We follow this approach here, 
as indeed we did in iDOd Dipole terms can be similarly treated, but are more important for the global fit studies. 
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Figure 3. Summary of template fit coefiicients. The top row presents results for the full-sky and individual hemispheres, the bottom 
row shows the values for the 33 regions (after excluding any pixels rejected by the EBV mask). For the regional plots, greater than 3-(t 
detections are shown with filled black circles whilst all other values are shown by open circles. The weighted averages at each frequency 
of the coefficients for the 33 regions are shown the filled red circles. Left column: Synchrotron - two fits for power-law spectral models 
are shown; in red the fit includes the 408 MHz reference point, in blue the fit only includes the 5 WMAP frequencies. Neither model 
includes sp ectral curvature. Middle co l umn: Free-free - the predicted coefficients (in intensity units) for different electron temperatures 
as given bv lDickinson. Davies fc David l|2003l ) are shown. Right column: Dust - Straight lines show power spectra fit to the data: in black 
is a single power law fit to K through Q; in red, two power laws are fit, one with a fixed index of /3 = 1.55. 



4 BASIC RESULTS 

At each WMAP band the emissivity of each of the 3 foreground components (free-free, dust and synchrotron) has been 
estimated as a ratio of template brightness; /xK R~^, ^iK/^Kfdss and /xK respectively. The analysis was a joint solution 
derived for all 3 components simultaneously, together with monopole and dipole terms. Results for the three templates can 
be found in Tables ETEl] and [C3l 

In Fig. 13] we provide a graphical summary of the results for both the global analyses and fits to individual regions. We 
present general inferences about the physical emission mechanisms below, and detailed comparisons of the fits to models in 
Section [S] 



4.1 Synchrotron 

As expected, the template fit coefiicients between the Haslam 408 MHz data and the WMAP sky maps all fall with frequency 
in a manner consistent with power-law emission. 

In Table [T] we define simple pairwise spectral indices between 408 MHz and the K-, Ka- and Q-bands. The global fits 
exhibi t a ty pical index of order —2.90 at K-band, with the results from all masks slightly fiatter than the value of —3.01 found 
in the iDOd analysis of the WMAP first -year data for t h e the n-pre ferred Kp2 sky- c overa ge. Analyses of the lower-frequency 
surveys at 408, 1420 and 2326 MHz by iGiardino et~aLl (|2002l ) and iPlatania et al.l l|2003l ) suggest a spectral index over this 
lower frequency range of approximately —2.7, thus our results support the idea of spectral steepening, continuing beyond 
K-band. We will consider this further in Section [S] 

However, there are differences in the coefiicients depending on the exact sky coverage that must refiect genuine spectral 
variations on the sky. The NFS is recognised to be an arc of steep spectrum emission at lower frequencies, thus it is not 
surprising that it is notably steeper than the rest of t he high latitude sky. The spectral index value of —3.03 at K-band is 
quite consistent with the value adopted in iFinkbeinerl ([2004) to remove the prominent emission from the WMAP data. The 
presence of the NFS also impacts measures of the spectral index in both the Galactic and Ecliptic northern sky, resulting in a 
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Table 1. Synchrotron fits between 408 MHz and WMAP K-, Ka- and Q-bands for various regions of the sky. For those regions that 
are detected at 2-cr significance, the fit values from Table ICll are converted into a spectral index using the usual power-law relation /3 
(Tf, CK u^). Associated errors are determined by using the fit values plus and minus the la error bars. Otherwise, one-sided 95% confidence 
level upper limits on the index are quoted based on the fit value plus 1.64cr. Bold text denotes detections at 3-(r significance or more. 
For the global fits, NPS - North Polar Spur, GN - Galactic North, GNjeduced - Galactic North with the NPS removed, EN - Ecliptic 
North, ENreduced - Ecliptic North with the NPS removed, GS - Galactic South, ES - Ecliptic South. 

modest steepening of the index. Interestingly, the northern Ecliptic hemisphere is notably steeper than the other hemispheres, 
and exhibits an increasingly steep index with frequency. Conversely, the corresponding southern hemisphere hints at spectral 
flattening, whereas both the north and south Galactic hemispheres are consistent with simple power-law behaviour. 

From the 33 regions of interest, there are 13 regions where the synchrotron fit coefficients are detected at greater than 3a 
confidence at K-band. Most of these regions include contributions close to the Galactic plane, although regions 2, 3 and 32 
are mostly at high Galactic latitude. The inferred spectral indices span the range —2.78 to —3.03, inconsistent with statistical 
variation alone, and likely representing genuine spectral variations on the sky. We note that region 3 contains the NPS and 
its behaviour seems to be dominated by that component. Significant emission is detected for region 8 at K-, Ka- and Q-band. 
Regions 13 and 14 correspond to an area of the sky containing the Gum nebula, and show no evidence for detection of 
synchrotron emission, with upper limits consistent with a steep spectrum, particularly in the southern region. Interestingly, 
region 9 has one of the steepest spectral indices on the sky, despite t he putative presence of the W MAP haze. However, the 
overlap between the brightest regions of the haze emission as seen in iDobler fc Finkbeineij (12008^ ) and regions 9, 16 and 17 
is very small, and unlikely to affect any studies here. 

It is interesting that the spectral indices inferred from the mean of the regional fit coefficients are steeper than the typical 
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global fit values. This trend was also seen in our previous work and may be due to a selection effect in that the 

regional subdivisions partly favour stronger synchrotron emission regions which may exhibit steeper spectra than normal due 
to synchrotron losses. In addition, evidence of spectral steepening with frequency is again generally seen although regions 10 
and 24 show inconsistent behaviour with the other regions in that a flattening of the spectrum is indicated. 

4.2 Free-free 

In this paper, we use an Ha template as a proxy for the free-free emission. We see significant correlation between the 
WMAP data and the template for the global fits at all frequencies (Table IC2|I . As with the synchrotron results, there 
are interesting variations depending on the exact sky coverage, with the northern Ecliptic hemisphere showing significantly 
enhanced amplitude, whilst the south indicates a lower emissivity. 14 individual regions are detected at Scr significance at 
K-b and. M ost of these regions lie close to the Galactic plane < 20°). 

IDDDI detail the relationship between the expected free-free brightness temperature and the related Ha intensity, and its 
dependence on both frequency and electron temperature (Te) in the ionised medium. As can be seen from their Fig. 5, the 
spectral dependence of the emission shows weak curvature, but over the range of frequencies covered by WMAP a reasonable 
approximation is a power-law with index —2.15. 

Inspection of the coef ficients in Table IC2I indicates that there are departures from this behaviour. Following 
iDobler fc Finkbeined l|2008bl ) we plot these results in intens ity units in Fig. [51 It should b e apparent that a bump in emission 
is seen around 40-50 GHz for the global fits, a feature that lDobler fc Finkbeineil (|2008bl ) argue is indicative of contributions 
from both classical free-free emission plus a spinning dust component in the WIM. The mean spectrum of the regional fits 
perhaps indicates a slightly broader bump in the spectrum. However, it is also the case that there is a range of behaviour seen 
amongst the individual regions, some of which are consistent with emission from a single physical component only - either 
free- free emission or a more steeply falling spectrum as expected from spinning dust. In Section [5l we will undertake a more 
detailed modelling of the emission in terms of these components. 

Such a feature in the spectrum of course has implications for the determination of physical parameters such as Te- 
Nevertheless, we can make some general inferences, particularly by examining the K- and W-band amplitudes that are least 
affected by a putative WIM spinning dust component. The global fits seem to be consistent with values in the range 6000 
- 8000 K, with some dependency on the exact mask used. The average free-free electron temperature inferred from the 33 
regions is also in this range, although for individual regions there is a spr e ad of values between 4000 and more than 10000 K. 
These values are somewhat higher than seen previously in iDavies et al.l l|2006l ). and this is due to the use of 3° smoothed 
data here, rather than the 1° resolution data used earlier for reasons provided in Appendix \K\ It is also inte r esting to note 
that the 7000 K temperature inferred from the 31.5 GHz CO-Bi?-DMR data at 7° resolution in lBandav et al.l (|2003i ) is quite 
consistent with the Ka-band values determined here. Moreover, the higher values in this paper are i n better agreement 
with the electron temp eratures deri ved from radio recombination line studies of extended HII regions (|Shaver et al.lll983l : 
iPaladini. Davies fc De Z otti 2004; A lves et al.| [2011') which derive an average Te value of ~ 7000 K in the vicinity of the solar 
neigh bourhood. There is therefore no need to invoke a large fraction of scattered Ha light to account for this discrepancy 
(e.g. , IWitt et aLllioiol ) . 

4.3 Impact of Dust Extinction 

iBennett etHI (|2003bl ') summarise various uncertainties in the use of the Ha template to trace free-free emission in the Galaxy. 
These include uncertainties in the value of the electron temperature Te, and in the value of the dust absorption correction, 
specifi ed here by the fd parameter. In principle, there will be variations in both throughout the Galaxy. As noted previously, 
Iddd ' have determined that for local regions such as Barnard's Arc and the Gum nebula, there is little absorption by dust, 
and this is the default situation t hat we have assumed in our analyses. However, for mid-Galactic latitudes a value fd ~ 0.3 is 
preferred, and iFinkbeinerl (|2003l ) adopted the assumption that the Ha emission is co-extensive with dust emission along the 
line-of-sight, ie. fd ~ 0.5. 

In this section we discuss the impact of varying fd on the template fit coefficients. In Fig. Uwe show this variation as 
a function of fd, and compare against the behaviour expected for a range of values of Te- In general, we expect that as the 
amplitude of the Ha template is corrected upwards by increasing fd, then the fit coefficients will decrease, implying a lower 
value for the electron temperature. This is indeed what is seen, but the extent of the correction depends on the region of sky 
under consideration. 

For the largest sky area that we consider, the EBV mask, the use of a template corrected absorption specified by fd = 0.5 
results in coefficients consistent with values of Te of order 3000 K. However, assuming no dust absorption correction yields 
values closer to 6000 K. A similar result is seen for the southern extension of the Gum nebula (region 14). If it is indeed the 
case that there is no evidence for dust absorption, then a higher temperature of 6000 K is determined. Conversely, regions 15, 
18 and 21 yield coefficients at K-band consistent with temperatures closer to 10000 K if no absorption correction is applied, 
whereas values of order 8000 K are found for fd = 0.5. Given that these regions exhibit rising spectra that may require the 
presence of significant emission from a spinning dust component, then the latter value would be more consistent. However, 
there are also regions of the sky typified by region 23, where apparently acceptable values of the electron temperature are 
inferred over a range of values of fd- Indeed, it may be that the large spread in coefficients seen for the different regions refiects 
changes in the fraction of dust mixed with the WIM as much as variations in Te- 
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Figure 4. Dependence of the correlation coefficients between the WMAP data and the template as a function of the absorption 
correction f^. The union mask refers to the combination of all 33 small regions defined in this paper. Thus the results for "Union + 
Galactic North" refers to the inverse-noise-weighted average coefficients for the subset of regions in the northern Galactic hemisphere. 
The observed dependencies for regions f3, f4 and 23 are representative of that seen for all regions with a significant detection amplitude 
at K-band. 



It does appear, therefore, that there are a range of values for Te and fd throughout the Galaxy, and reaching conclusions 
about their values solely from studies of He correlations suffers from degeneracies be t ween t he parameters. The adoption of a 
value for the dust absorption correction of 0.5 in, for example. iDobler fc FinkbeineJ (|2008bh . thus seems to be associated, at 
least in part, with the low electron temperatures inferred. This uncertainty then has implications for modelling the emission 
of the diffuse component. 

A more serious compli c ation would arise if the validity of using Ha as a tracer of free-free emission were questioned. 
iMattila. Juvela fc LehtinenI (|2007l ) have argued that the H^ excess towards the high-latitude interstellar cloud LDN 1780 
is the result of scattering of H^ photons produced elsewhere in the Galaxy by dust in the cloud. Earlier computations by 
IWood fc Reynolds! l| 19991 ) suggested that this contribution would typically be 5-20% at high-la titudes, alth ough their model 
has been criticised due to the assumption of a smooth distribution of material in the ISM. et al. 1 |2010i l have used 

an empirical relation to relate the scattered Hq intensity in the translucent cloud LDN1780 to the IRAS 100 fim diffuse 
background intensity and conclude that this estimate is reasonable for 50% of the high-latitude sky, but that the scattered 
contribution can be highly structured and result in contributions of between 25 and 50% of the observed intensity for a further 
40% of this region of the sky. Such a result would clearly have implications for using an He template to trace free-free emission, 
and a complex relationship between the template and dust would arise both due to the scattering contribution and due to any 
dust absorption correction applied to the data. Applying a correction for the former would effectively result in an increase of 
the Hq- WMAP correlation coefficients and the refore of the i nferred electron temperature, whereas application of the latter 
results in the opposite behaviour. More recently. iB randt fc Dr ainc (20 1l|) have estimated th e fraction of high-latitude H^ that 
is scattered to be 19 ± 4%, a value consistent with that proposed bv lDong fc Drain3 (|201ll ) to reconcile the low ratio of radio 
free- free to H^. 

Ultimately, unravelling the degeneracy between the electron temperature. He scattering and dust absorption requires 
additional observations. Detailed RRL su rveys in th e Galactic plane together with radio-continuum surveys at frequencies of 
~5 GHz, as expected from the G-BASS (.King et al.i , 20101 will be important in this respect. 

Finally, we would like to make some remarks about region 11, the coefficients of which suggest an exceptionally high 
temperature of more than 25000 K. If this result were considered unphysical, then naively a value of fd ~ 1 would be required 
in order to lower the inferred temperature to the 6000 K seen in the EBV fit. In fact, a more realistic assessment of the 
situation is that the EBV mask is not large enough to eliminate some parts of region 11 close to the Galactic plane where the 
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simple dust absorption correction is untrustworthy. If instead we apply the KQ85 mask before analysing the region, then for 
a more plausible value of fd = 0.5 the inferred electron temperature is again consistent with 6000 K. 

4.4 Dust 

As can be seen in the right-hand panels of Fig. (3] the template fit coefficients determined between the WMAP data and the 
FDS8 dust template prediction at W-band are consistent with a rising thermal dust contribution at frequencies higher than 
V-band, and a rising spectrum to lower frequencies below it. The latter corresponds to the now widely identified anomalous 
microwave emission (AME). 

It appears that the FD S8 te mplate underpredicts the W-band amplitude by approximately 30% for the global fits, 
consistent with the results of iDOfil . A broad range of values for the individual regions is seen, but only three detect emission 
at a statistically significant level. The mean emissivity at W-band of all regions shows a more modest enhancement, but is 
nevertheless consistent with the FDS8 predictions, as indeed are the three significant regions. 

At K-band, there are variations in the global fit amplitudes depending on sky coverage. The KQ85 mask indicates a lower 
emissivity compared to the EBV as might be expected. The Southern Galactic and Ecli ptic h emispheres have the highest 
amplitudes. All coefficients are higher than those determined for the Kp2 sky coverage in 'DOfj. 26 of the individual regions 
detect emission at 3(t sign ifican ce or higher. The mean emission amplitude lies between that for EBV and KQ85, and again 
somewhat higher than in The regions indicate a variation around the mean of approximately 50% of its amplitude, 

inconsistent with statistical errors alone and indicating genuine spatial variations in the AME emissivity. 

The overall spectrum would appear to be well-described by a superposition of two power-law emissivities. The thermal 
dust emission described by the FDS8 model is adequately represented by an emissivity index Pd ~ 1-55 over the WMAP range 
of frequencies. We fit the power-law AME spectrum with the thermal dust index fixed to this value, and find AME spectral 
indi ces of order —2.7 fo r both the global fits and regional mean. 

iDraine fc LazarianI (|l998al ) first proposed that the AME could be explained by electric dipole radiation from rotationally 
excited small interstellar grains, or spinning dust. Our results have ramifications for such models of the emission. In particular, 
given that the spinning dust spectra typically fall off steeply with frequency beyond their peak, then it is unlikely that a single 
such spectrum could account adequately for the effective power-law emission. Indeed, the observed spectrum is presumably 
formed from a superposition of components with varying spectra as a consequence of their differing physical environments 
along a given line-of-sight. We will discuss detailed fits of the observed emission in Section [SI and their implications for more 
refined models of the AME. 



5 MODEL FITS 

In order to compare the derived template fit coefficients of a given foreground comp onent with various theoretical ly motivated 
foreground models, we adopt the simple procedure used previously in, for example. iDobler fc FinkbeineJ l|2008bl ). The model 
parameters of the foreground model are extracted using the least-square minimization defined as, 

where F{vi) and (j{vi) are the observed fit coefficients and standard deviation at WMAP frequency bands and F^' [v) axe 
the fit coefficients given the foreground model. We use a Levenberg-Marquardt method to determine the coefficients, and 
quote the reduced chi-square for a given model fit. The degree of freedom (dof) is defined as n — m where n represent the 
number of data points to be fitted and m represents the number of free parameters for a given foreground model. We note 
that the interpretation of such values can be problematic when the number of data points is small, as is the case here. This is 
compounded by the correlated errorbars that are dominated by a CMB common to all frequencies. Indeed, the model fits are 
only to be considered as indicative rather than definitive. A more robust approach would apply a multi-frequency analysis to 
account for the correlated errors at each frequency due to the dominant CMB term in the covariance matrix. This will form 
the basis of a future publication where models are fitted directly to the data rather than to previously derived template fit 
coefficients. 

5.1 Spinning Dust preamble 

lEricksonI (|l957l ') first pro posed the basic mechanism of spinning dust emission from the rotation of small dust grains with 
electric dipole moments. iFerrara fc Dettmaij l|l994l ) later suggested that such grains in the diffuse ioni sed medium should 
exhibit significant radio emission peaking at a frequency between 10 and 100 GHz. However, it was iDraine fc LazarianI 
lll998al) that sugge s ted suc h emission could explain the AME, and provided detailed computations of its spectral shape in 
IDraine fc LazarianI (|l998bh . 

Indeed, testing for the presence of such emission requires detailed pre dictions of the emission spectra to com pare with 
observations. However, given the large number of parameters in the model l|Ali-Haimoud. Hirata fc Dickinsonll2009l ). and the 
fact that most lines-of-sight likely average over many emission regions, this is difficult unless specific objects are considered. 
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Such an analysis was performed r ecently with early data from the P lanck mission iTauber et al.ll2010l ). specifically for the 
Perseus and p-Ophiuchus regions jPlanck Collaboration et al.ll2011dh . Thus, it is usual to adopt spectra for spinning dust in 
a variety o f phases of the ISM, computed for 'typical' values of the parameters for those physical conditions. Alternatively, 
iGold et al.l (|200^) proposed a generalis ed analytic form for the spinning dust emission for their analysis of the WMAP 5- 
year data (see also iBonaldi et al]|2007l ). However, a fit to the exact CNM form from DL 98 underestimates the emis s ion at 
frequencies beyond the peak. This is problematic given that more recent calculations bv lHoang. Draine fc LazarianI (|2010l ) 

indicate that the emission in this region may be broadened by the inclusion of addi tional physical processes. 

We will adopt spinning dust emission templates based on the SPDUST2 code (jSilsbee, Ali-Hai'moud fc Hiratalboi l) for 
emission in the cold neutral medium (CNM), warm neutral medium (WNM), and warm ionised medium (WIM). We then 
allow the spectra to be shifted in both amplitude and peak frequency to fit the data, considering that this approximately 
mimics the effect of varying the physical parameters as actually required. 



5.2 Synchrotron 

The following analytic forms are used to fit the synchrotron coefficients. 

• Model SI : Given the power law energy distribution of cosmic ray electrons, we assumed a power-law emissivity in terms 
of brightness temperature over the WMAP frequencies as 

Ta{v) = As >^ (Q^' (4) 

V23/ GHz 

where j3s is the spectral index and As is the normalised amplitude expressed in /iK with respect to the frequency vo = 23 GHz. 

• Model SII : A power-law emissivity is assumed to extend from 408 MHz up to and through the WMAP frequencies. 
Since we use the 408 MHz survey as a template for the synchrotron emission, the amplitude at the low frequency must be 
reproduced perfectly. This results in an effective constraint to be applied to the analytical form above, and we then fit for the 
spectral index /3s only. 

Ta{v) = 10'' (5) 
V 0.408 /GHz 

• Model SIII : The cosmic ray electron energy spectrum is expected to steepen with time due to ra d iation energy 
loss. A review of c osmic-ray propagation including electrons can be found in IStrong. Mo skalenko fc PtuskinI (|2007l ). whilst 
IStrong. Orlando~&: Jaffe (2011i ') directly test propagation models based on cosmic-ray and gamma-ray data against synchrotron 
data from 22 MHz to 94 GHz as averaged over mid-latitude regions (10° <| 6 |< 45°). The latter analysis confirms the need 
for a low-energy break in the cosmic-ray electron injection spectrum to account for the steepening synchrotron spectrum. 
Since we do not include s ynchr otron information at frequencies intermediate to 408 MHz and the WMAP data, we follow the 
treatment of lCold et al.l ( 20091 ) . Specifically, the emissivity is assumed to follow a power- law from 408 MHz until K-band and 



then to exhibit spectral curvature as follows, 

TAiiy) = As X {-^^ V<VK 



' GHz 



V 



= ^» ^ ( ^ ) V>PK 
V 23 /GHz 

As above, the 408 MHz point is fixed, thus As can be written in terms of /3s, and we are left to fit this spectral index and the 
curvature Cs. For a WMAP frequency point, the above equation reduces to the form: 

T.(^) = 10«x(-^)^^ x©"""^' (6) 
V 0.408 /GHz V 23 /GHz 



The results are summarised in Table [21 Fig. [5] presents a comparison of the model fits for three of the global masks, 
plus three of the regions that can be considered representative of the general results seen. The results generally should be 
considered indicative rather than definitive, as would be expected given the small number of data points and the comparable 
number of model parameters. 

The synchrotron fit coefficients obtained for the global masks that include the 408 MHz datum as a reference point 
(model SII) indicate a typical power-law spectral index of Ps = —2.9 ± 0.1. Model SI fits to the 5 WMAP frequencies are 
generally steeper, however, only the Ecliptic North region shows evidence of such behaviour at a significant level (a 2(t shift 
in the spectral index). Indeed, a model with negative spectral curvature is the best fit to the data, indicating steepening due 
to cosmic ray energy loss mechanisms. Curiously, the Ecliptic South indicates spectral fiattening, although at much lower 
significance. The SII fit to the NPS region finds a sp ectral index Ps = —3 .05 ± 0.03 which is steeper than the average spectral 
index variation over the remaining sky, as expected (|Lawson et al.|[^987^ . 

For the individual regions, the results are generally similar to the global fits. The spectral behaviour is consistent with 
emission from power-law cosmic ray electron spectra with only hints of a steepening between 408 MHz and K-band. Region 11 
might be considered to show weak evidence for negative spectral curvature. More interestingly, three regions (10, 12 and 24) 
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Table 2. Model fits to the synchrotron coefficients determined between the 5 WMAP frequencies and the Haslam 408 MHz template for 
large sky areas (upper part of table) and for those regions that indicate a Scr significant amplitude at K-band (lower part). The models 
SI, SII and SIII are described in Section 15.21 As represents the normalisation amplitude at K-band, /3s the synchrotron spectral index 
and Cs the spectral curvature. The key for the global fits is as for Table [T] 
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Figure 5. Synchrotron spectral fits for various regions of the sky. A comparison is made of the observed template fits amplitudes derived 
from the 5 WMAP frequency bands and the 408 MHz survey with the 3 models SI, SII and SIII as defined in Section 15.21 The observed 
spectral dependencies of regions 3, 10 and 15 are representative of those seen for all regions with a significant detection at K-band. 
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indicate positive curvature that might be understood as due to the presence of multiple emission regions with varying spectral 
behaviour along the line-of-sight. However, region 10 has a very fiat spectrum between K- and W-band that is strongly 
favoured over a steeper spectrum extending to 408 MIfz, and perhaps even mildly inconsistent with the curvature model 
adopted here. Whether this indicates some problem with the analysis in this region, due to problems with the templates or 
cross-talk between components, is difficult to determine. It is apparent that the region lies close to the North Celestial Pole 
where the 408 MHz template still retains considerable striations from the original survey, whereas we certainly consider that 
the template fitting methodol ogy has by now be en extensively tested. However, the spectral flattening seen in other analyses 
is less dramatic. For example, I Gold et al.l ('2011) find a range of values of Cs between ±0.7 with a variance of ~ 0.4 for their 
MCMC analysis of the WMAP data combined with the 408 MHz data. 

5.3 Free- free 

The free-free spectral index is almost constant over the WMAP frequency range and changes only slightly with the the electron 
temperature. Over this frequency range, the spectral dependence of the corresponding Hg emi ssion can then be approximated 
with a power-law model of fixed spectral index —0.15. However. iDobler fc Finkbeineij (|2008al ) observed deviations from such 
power-law behaviour, which they attributed to the presence of spinning dust emission in the WIM also traced by Hq. We 
investigate our template fit coefficients in terms of the following three models. 

• Model FI : We consider that the emission is due entirely to the free-free mechanism, ie. a power-law model emissivity 
with a fixed spectral index of —0.15 is assumed. 



Thus only the free-free amplitude [Af) needs to be estimated from the template fit coefficients. 

• Model FH : The coefficients are fitted with an empirically motivated power-law model emissivity. This is particularly 
illustrative as to the extent that the standard free-free emissivity index is a poor fit to the data. 



In this case, both the foreground amplitude {Ap^) and spectral index [Ppl) are to be estimated. 

• Model Fill : The coefficients are fitted with a comb ination of free-free emission and a spinning dust model computed by 
the SPDUST2 code (|Silsbee. Ali-Haimoud fc HiratalbOllT ) for typical WIM conditions. This model peaks around 30 GHz, but 
we allow a simple shift Azavim to be applied to the spectrum, as an approximation that represents the effect of varying the 
WIM physical parameters to match the model spectrum to the data. 



Here, DwiM(i^) represents the normalised spinning dust spectral model for the WIM at a given frequency. Clearly, we 
must now fit 3 parameters - the free-free amplitude {Af), WIM amplitude (Awim ) and WIM frequency shift (Ajavim). 
We have also considered fits of this model to the coefficients derived using the He after correction for dust absorption [fd ~ 0.5). 



The results are summarised in Table O Fig. [H] presents a comparison of the model fits for three of the global masks, plus 
three of the regions that can be considered representative of the general results seen. 

The global fits all indicate a deviation from the free-free emission law (FI) at high significance. Unconstrained power-law 
fits (FII) are generally fiatter, and models containing both free-free and spinning dust emission (FH) are typically preferred. 
For the EBV mask, the ratio of free-free to spinning dust emission (at the shifted peak frequency) is of order 7, but this drops 
to 4 for the Galactic North. All regional fits are consistent with a 15 GHz shift to higher frequencies of the peak intensity for 
the WIM spinning dust component with the exception of the Galactic North region which favours a smaller value of 5 GHz. 

The fits to individual regions are generally consistent with the global mask results, and only region 9 shows a clear 
rejection of model FI. The typical frequency shift is again of order 15 GHz, with a larger dispersion with several regions 
preferring no shift at all. Of those, regions 7 and 12 can be explained by a single component - either free-free, power-law 
emission or spinning dust alone with no frequency shift. However, it should be noted that the H^ fit coefficients for these 
regions are only significant at K-band. Region 11 is consistent with free- free only emission, and yields no evidence for spinning 
dust, but the region itself has been flagged as anomalous as discussed previously in section 331 

Region 13 corresponds to the Northern Gum nebula, and is dominated by free-free emission with a small contribution 
from spinning dust at a lower peak frequency, ^ 35 GHz, than is typical. The Southern part of the nebula is contained in 
region 14. This is the most clearly detected structure on the sky, significant at all frequencies, as traced by the Hq. template, 
and shows significant evidence for a spinning dust contribution with a free-free to spinning dust ratio of ~ 5. Otherwise, this 
ratio varies considerably from region to region. 
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Table 3. Model fits to the free-free coefficients determined between the 5 WMAP frequencies and the DDD Hq template for large sky 
areas (upper part of table) and for those regions that indicate a 3(t significant amplitude at K-band (lower part). The models FI, FII and 
Fill are fully defined in Section 15.31 Af represents the normalisation amplitude of the free-free emission at K-band, Api^ represents the 
normalisation amplitude for the power-law model emission at K-band, f}pL the corresponding power-law spectral index, AwiM is the 
amplitude of the WIM spinning dust model in units of lO'^" R cm~^, and At^jM the shift in frequency of the peak of the dust model 
to better the fit the data. The key for the global fits as for Table [T] 
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Figure 6. Free-free spectral fits for various regions of the sky. A comparison is made of the observed template fits amplitudes derived 
from the 5 WMAP frequency bands and the DDD H^j template with the 3 models FI, FII and Fill as defined in Section 15.31 Also 
shown in the plot are the separate free-free (FF) and WIM spinning dust components that constitute model Fill. The observed spectral 
dependencies of regions 7, 13 and 14 are representative of those seen for all regions with a significant detection at K-band. 
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Region 
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5100 ± 2200 
6000 ± 200 
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3700 ± 600 
5900 ± 2900 
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2300 ± 900 
3200 ± 300 
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3300 ± 1300 
3500 ± 1500 
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7700 ± 1700 
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4000 ± 1800 
5400 ± 1000 
2300 ± 700 
5600 ± 2500 
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Table 4. Inferred free-free electron temperature Te in Kelvins corresponding to model fit Fill with dust corrections fd = 0.0 and 0.5. 
(a) These regions show an effective degeneracy between free-free only or WIM spinning dust only solutions, thus Fill solutions have no 
contribution from free-free emission. We have therefore used the FI results to compute Tg in these cases. 

Some regions that have a rising spectrum in terms of the coefficients and naturally favour models with a spinning 
dust contribution over pure free- free are better fitted still by a flatter/rising power-law emission model. This might be 
alleviated with more detailed spinning dust modelling (beyond the scope of our paper), or by including physical effects that 
increase the spinning dust ampl itude at frequencies higher than the peak. Such modifications have been investigated by 
iHoang. Draine fc L azarian (2010'). 

The A/ fit coefficients can, of course, be converted to estimates of the thermal electron temperature in the ionised 
medium. Since the global mask fits require the presence of a spinning dust component, we only consider the Te results from 
model FIJI, as presented in Table 21 The table also includes results derived from fits to an template corrected for dust 
absorption assuming fd = 0.5. We do not include the detailed coefficient results here since the interpretation presented above 
remains essentially unchanged, and only conclusions about Te are affected. The global masks are consistent with values of 
the electron temperature of order 6000 K without any dust absorption correction, falling to 3000 K when fd = 0.5 is assumed. 
These values are for guidance only - since the regions are not independent, an average is meaningless. It is interesting to note 
that the Ecliptic North shows an enhanced temperature some 50% higher than these typical values. Whether this reflects 
some property of the Ha template is difficult to say, but the coverage is dominated by measurements from WHAM data. 
Conversely, the Ecliptic South global fit gives a lower temperature that may also suggest issues with the template, given 
that it is largely co mprised of the SHASSA fields. New observations from the southern extension of the WHAM survey 
iHaffner et al.|[2010l l should help resolve this issue in the near future. Nevertheless, problems may still remain near the ecliptic 
poles given difficulties with removing the geocoronal Ha contribution. For the individual regions (ignoring region 11 which 
is considered to be anomalous), we find a weighted average of 6300 ± 200 K (dropping to 5900 K if region 11 is excluded) 
without a correction for dust absorption, and 2900 ± 100 K otherwise. These values are in good agreement with the global 
averages as might be expected. However, note that the dispersion of values is considerably larger than the quoted uncertainty, 
implying true variations in temperature on the sky. Moreover, as discussed previously, it is likely that some of the dispersion 
seen refiects the existence of a range of values for both Te and fd throughout the Galaxy. 



The total dust emission is modelled as a combination of the relatively well-understood thermal dust emission and the AME. 
The former is assumed to have a fixed spectral index relative to the FDS8 94 GHz template over the WMAP frequency range 
as determined directly from the FDS8 model. We consider the following two models in order to fit the dust coefficients. 

• Model DI : The dust coefficients are fitted with a combination of thermal (vibrational) dust and a power law dust- 
correlated AME. 



Three parameters - the thermal dust amplitude {Atd), the power law dust amplitude {Apld) and power law dust spectral 
index [Ppld) - are fitted to the coefficients. 



5.4 Dust 
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Table 5. Model fits to the dust coefiicients determined between the 5 WMAP frequencies and the FDS8 template for large sky areas 
(upper part of table) and for those regions that indicate a 3(t significant amplitude at K-band (lower part). The models DI and DII 
are fully defined in Section 15.41 Api^]j and Ppld represent the normalisation amplitude at K-band and spectral index respectively of a 
power-law anomalous dust emission component. Ay^NM is the amplitude of the WNM spinning dust model at K-band, Aqnm is the 
amplitude of the CNM spinning dust model normalised at 41 GHz, and AvcNM is the shift in frequency of the peak of the CNM dust 
model to better the fit the data. Atd is the amplitude of the thermal dust emission with an assumed spectral index Ptd = 1-55 as 
determined directly from the FDS8 dust model. The key for the global fits as for Table [Tl 



• Model DII : The dust coefficients are fitted with a combination of thermal dust and two spinning dust components 
( "CNM and WNM) The t wo spectra are generated using the SPDUST2 code assuming typical CNM and WNM conditions 
iDraine fc Lazarianlll998"al ). Given that both spectra peak at approximately the same frequency (« 30 GHz), and the limited 
number of degrees of fr eedom available in the fit, i t is on ly possible to apply a frequency shift to one component in order 
to match observations. iHoang. Lazarian fc Draind l|201ll ) found that modifying the CNM properties to increase its peak 
frequency yields a closer match to the WMAP observations, and therefore we elect to allow a frequency shift of this component. 

TAiy) = ^WNM X DwNm(!^) + ^CNM X Dc^yy(y - A/^CNm) + AtD X — (11) 

V 94/ GHz 

We fit four parameters to the template fit coefficients: the thermal dust amplitude (Atd), the WNM amplitude (Avvnm) 
normalised at 23 GHz, the CNM amplitude (Acnm) normalised at 41 GHz, and the CNM peak frequency shift (Ai/cNM). 



The results are summarised in Table [5] Fig. [7] presents a comparison of the model fits for three of the global masks, plus 
three of the regions that can be considered representative of the general results seen. 

It should be apparent from the table that both models DI and DII provide adequate fits to the global masks. However, in 
order to explain the emission at V-band in terms of spinning dust models, two contributions are required in combination with 
an enhanced thermal dust normalisation some 20% larger than that predicted at 94 GHz by FDS8. Such models also require 
a shift in the peak CNM emission of ~22GHz and a ratio of WNM to CNM spinning dust emission of order 4:1. In fact, if 
the CNM emission were to be omitted, or the thermal dust amplitude kept at the canonical FDS8 value, then the V- and/or 
W-band amplitudes would be substantially underpredicted. In practise, given the low number of data points sensitive to the 
thermal dust emission, there is an effective degeneracy between the CNM and vibrational dust components, that may be 
alleviated somewhat if the spinning dust models can be enhanced in amplitude at frequencies higher than their peak. Indeed, 
spinning dust can contribute significantly to emission in the ^^60 to 90 GHz range depending on the local physical conditions. 
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Figure 7. Dust spectral fits for various regions of the slsy. A comparison is made of the observed template fit amplitudes derived from 
the 5 WMAP frequency bands and the FDS8 template with the 2 models Dl and DII as defined in Section 15.41 Also shown in the plot 
are the separate power-law anomalous dust (PLD) and vibrational (thermal) dust (VD) that constitute model Dl, and the cold neutral 
medium (CNM) and warm neutral medium (WNM) components that, together with the vibrational dust, constitute model DII. The 
observed spectral dependencies of regions 12, 14 and 21 are representative of those seen for all regions with a significant detection at 
K-band. 



However, one might then expect to see a much flatter spectrum around the Q-, V-, and W-bands than is actuaUy observed. 
It may also be that some modification of the assumed thermal dust spectral index is required. 

The fits to individual regions are less informative, as expected given the lack of clear detections for the template fit 
coefficients at V- and W-band. Given this, the analysis might be criticised for overfitting the data points. Nevertheless, both 
a power-law and a WNM spinning dust component are good candidates to explain the lower frequency points. In almost all 
cases, a significant contribution from CNM or thermal dust emission is not required. This discrepancy between the regional fits 
and global analysis may simply refiect a signal-to-noise issue at higher frequencies than the former, or the fact that the latter 
are essentially averages over many regions which exhibit some variation in dust properties. Regions 12 and 16 are interesting 
in that the power-law model is clearly preferred over model DII, although the analysis is really only constrained by the data 
points at K- and Ka-band. For those regions in which a significant vibrational dust contribution is required, the amplitude is 
systematically higher than the canonical FDS8 value, although still consistent with it. However, the required Dl enhancement 
is lower, thus the DII results more likely again point to a deficiency in the specific spinning dust models used in this analysis. 

The FDS8 model value for the dust spectral index used here is quite fiat, whereas early Planck papers 
l|Planck Collaboration et al.l 2011bl al) studyin g emission at frequencies higher than ~100 GHz favour values closer to 1.8. 
However. iPlanck Collaboration et al.l (|201ld) suggests that the dust SED flattens in the millimetre wavelength range. At 
lower frequencies, there is considerable uncertainty, and detailed modelling for candidate AME regions (including Perseus and 
p-Ophiucus) of the foreground spec tra for all components and across the Planck frequencies suggests values in the range 1.5-2 
(|Planck Collaboration et al.ll2011dl ). We have repeated the calculations above but imposing a thermal dust spectral index of 
1.7 on the fits. For model Dl, there is a general increase in amplitude of the power-law dust component, with an associated 
flattening of the spectral index at levels below the 1-a error bar, and no impact on the thermal dust amplitude. Model DII 
shows similarly modest increases in the amplitudes of the two spinning dust components, with little change in the CNM peak 
frequency shift, and again no change in the thermal dust amplitude. However, adopting the FDS8 thermal dust model at 
W-band and then imposing a steeper index to lower frequencies is not self-consistent. Therefore, a more detailed treatment 
is required including specific modelling of all dust components combined with higher frequency measurements. An obvious 
potential error in our analysis is the selective use of the CNM and WNM spectra for typical conditions in those phases of the 
ISM. What is clear is that the frequency range of 60-100 GHz is a remarkably interesting regime for dust astrophysics. 
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Table 6. R.m.s amplitudes of the foreground templates. These values are scaled to the WMAP frequencies using the fit coefficients in 
Tables [CTI [C2l and fel 



6 DISCUSSION 

In our study of the free-free, dust and synchrotron foreground components in the WMAP data, we have selected fields based 
on the morphological properties of three templates that trace the emission at wavelengths where the emission mechanism is 
(largely) dominant, and we assume that such partitions correlate well with spectral behaviour. Each of these components has 
then been quantified in terms of a mean value of the emissivity in each of the 5 WMAP bands. 

In Fig. [S] we show the r.m.s contributions of the foregrounds in antenna temperature as traced by the templates for three 
of the global fits and a representative set of the regional results. These values were derived from the template r.m.s values as 
recorded in Table[B]and scaled appropriately by the fit coefficients from our analysis. For comparison, we show the uncorrected 
r.m.s. amplitude at each WMAP frequency, together with an estimate of the CMB fluctuations from the ILC sky map, all 
determined for the same sky coverage as the foregrounds. We assume that this template based technique tr aces essentially al l 
of the foreground contributi on, and this may not be the case. In particular, we recognise the WMAP-haze l|Finkbeinerll200j : 
iDobler fc Finkbeinedr2008ah as an important exception. Nevertheless, our interpretations should remain robust. 

Note that the spectral shape of the integrated foregrounds, where determined, is relatively simple when sparsely sampled 
by the WMAP frequencies, but nevertheless it is interesting that a simple power-law frequency dependence, or only modest 
devia tions therefrom , is n eeded to describe the emission from K- to V-band. This observation was central to the work of 
iPark. Park fc GottI | 200?! ) in partitioning the sky to perform local CMB reconstruction. 

For the EBV mask, it is clear that the foregrounds as traced by three templates are of similar amplitude at K-band, 
and the free-free and dust contributions remain roughly equal to V-band, whereas the synchrotron falls off more steeply. The 
foregrounds remain comparable to, or larger than, the CMB fluctuations to Q-band, and significant at V- and W-bands. 

Comparing the total foreground contributions between the North and South Ecliptic hemispheres, it can be observed 
that the amplitude is higher in the South at least from K- to V-band. It also appears that the free-free component shows a 
particular enhancement in the South. Whether this may be connected to problems in the template (see Appendix \K\ for 
details) might be discussed. However, the foreground differences are interesting in the conte xt of the hemispherical asymmetry 
seen in the distribution of power in the CMB in Ecliptic coordinates, as noted originally bv lParkI (|2004 ): lEriksen et all (|2004|[ ) 
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and revisited many times in the literature (see ICopi et al.|[2010l . for a review). Indeed, the r.m.s. CMB signal plotted here 
clearly shows this asymmetry. Specifically, if the foregrounds were underestimated in the Northern hemisphere (as might be 
the case for the free-free component) , then correcting for this effect would only exacerbate the problem. A similar consequence 
results from an overestimate of the foregrounds in the South. Thus, it seems unlik e ly tha t the CMB asymmetry is connected 
with problems of Galactic foreground estimation, as noted before in iHansen" et al.1 l|2006l ). 

As expected, there are also clear variations between the regions used in our study. Region 9 represents a region with 
strong foreground contributions at all frequencies, the dominant contribution being from the free-free emission. Region 23 is a 
dust dominated region, but the significance of the contamination relative to the CMB is lower, although the dust is detected 
clearly at all frequencies. Finally, region 33 typifies many high- and mid-latitude regions where foregrounds are only detected, 
if at all, at K- and possibly Ka-band. The amplitude of the 95% upper limits combined with the agreement between the 
uncorrected frequency maps and the CMB estimate from the ILC implies that foregrounds are essentially unimportant in 
these regions. 

There are several interesting aspects of the study with regards to the properties of the foregrounds that should be 
highlighted. The synchrotron emission shows clear evidence of steepening relative to measures at GHz frequencies. There is a 
hint of this steepening continuing beyond K-band. The temperature of the electrons responsible for the free-free emission is 
now much more consistent with other indicators such as RRLs, ie. ~7000 K at high latitudes. This result is in part due to an 
improved understanding of the Ifa template. This clearly needs to be improved and will benefit greatly from the WHAM-S 
survey when available. However, uncertainties related to the dust absorption remain. The AME is essentially ubiquitous, 
comparable in amplitude to the free-free emission, and with a spectrum broader than a single component. There does appear 
to be an additional anomalous component correlated with the H^ emission, but there are variations in relative strength on 
the sky. An improved model of the AME requires an improvement in the understanding of the corresponding thermal dust 
contribution. Specifically, the emissivity of the thermal dust tail to low frequencies must be determined, and a detailed model 
of the emissivity as a function of wavelength and dust temperature needs to be derived from the data. 

Finally, we note that studies of foregrounds at microwave wavelengths using the techniques in this paper are essentially 
limited by the quality and resolution of the available templates. A comparison of synchrotron and free-free properties over the 
full-sky are, at best, possible on angular scales of 1° and above, given the resolutions of the 408 MHz and H^ data. Moreover, 
we have noted some problems with the H^ template on scales less than 3°. However, this should not limit the analysis 
of Planck data, and especially that from the high-frequency instrument (HFI) operating at frequencies above 100 GHz, 
which will determine the R- J tail of thermal dust emission. Similarly, data at frequencies ~ 5 — 20 GHz will help determine the 
complicated mix of low frequency foregrounds. For example, there could be a possible contribution from synchrotron radiation 
with a fl atter s pect ral index compared to the values determined from correlations with 408 MHz data. We note that a recent 
study bv lPeell ([2011) who used a 2.3 GHz survey of the southern sky found that the AME was relatively unaffected limiting 
the contribution of such a component to ~ 7% of the AME at 23 GHz. 
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Figure 8. R.M.S amplitude in antenna temperature observed for various regions of the sky. The signal associated with the total 
foreground is shown as black filled circles. The synchrotron contribution traced by the Haslam template is represented by red squares, 
the free- free traced by the Hq template as dark blue triangles, and the dust contribution traced by the FDS8 template as purple diamonds. 
One-sided 95% c.l. upper limits are shown for non-detections, with the appropriate colour coding for the component represented. The 
r.m.s. amplitude of the uncorrected WMAP frequency maps is indicated by open light blue circles. The orange line corresponds to the 
CMB level estimated from the 7-year ILC map. The observed spectral shapes of regions 9, 23 and 33 are representative of those seen for 
all regions. 
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Figure Al. Left and middle: difference map between tlie IDDdI and IFO j Hq templates. Tlie iFoj map lias been additionally smoothed 
by a Gaussian beam of (60^ — Q'^)^/'^ degrees so that the two sky maps have a matched 1° resolution. The two plots are identical except 
for the color scale. Right: the coverage of the different datasets used in the two te mpla tes. In light green is where both templates use 
WHAM dat a, and in red where both use S MAS S A data. Orange shows where the iFOa template uses VTSS data. In dark blue is the 
region where IDDdI uses WHAM data while IFOSI uses SHASSA. 
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Figure A2. Pseudo-C^ spectra computed from the IDDdI and IFO j Ha templates. Left: pseudo-C^ spectra of discs centered on either 
Celestial North or South (|(5| > 2 0°). Dashed lines show the Southern spectra while solid the Northern. Results from the IDDO template 
are shown in red, those from the lFOa map smoothed to 1° resolution are in black, and in orange we show results from the lFOa template 
without additional smoothing. Dotted lines give simple power law spectra smoothed with Gaussian beams of 1° (1.3°) in blue (green) 
with arbit rar y powe r law indices and amplitudes fit by eye solely for visual comparison. Middle: pseudo-Cf spectra of sky regions where 
both lFOsI and IDDdI templat es us e only WHAM data. The | F03| template has been smoothed as described in the text. Right: pseudo-Cf 
spectra of sky regions where IFOSI uses VTSS data but lDDpf uses WHAM. For comparison, the pseudo-C^s for the masks themselves are 
shown, the VTSS used for the red and black curves, and the North mask used in the top figure. The effects of the highly structured 
VTSS mask dominates toward high £, which makes the lower two plots difficult to compare with the top plot, where an azimuthally 
symmetric mask does not impact the power spectrum. But it is apparent that the differences in the north are due to the use of the VTSS 
data. 



APPENDIX A: H„ TEMPLATES 



We have compared the IDDDI and iFinkbeinerl (|2003l . hereafter F03) Hq, templates and performed template fits to the 
WMAP data using both for comparison. As in iDOd there are sometimes significant differences in the resulting cross-correlation 
coefficients, that we previously ascribed to an inaccuracy in our knowledge of the WHAM effective beam. Here, we investigate 
the detai led dif ferences between the templates themselves. 

The IDDDI template used a combination of W HAM data in the north and SHASSA data in the south, with preference 
for the WHAM data down to a declination of — 10°. If03I chose to u se SHASSA data wherever i t was available, augmented by 
data from the Virginia Tech Spectral-line survey (hereafter VTSS, iDenn ison. Simonetti fc Tbp asna 1998), a complementary 
narrow-band imaging survey which covers the northern sky for 5 ^ —15°, and has an angular resolution of 1.6 arc minute. The 
SHASSA data wer e proc essed differently, and different methods used to adjust the resol ution s of these datasets for consistency 
and merging. The IDDdI template claims a nominal resolution of 1° FWHM, while the lF03l template has a stated resolution 
of 6'. 
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Fig. lAll shows the difTerences between the two maps on two different color scales to bring out the small and large amplitude 
differences. Clearly, the processing differences result in maps that differ at the roughly 1 Rayleigh (R) level over large regions 
of the sky, and of more than 20 R level near very bright regions, although the latter are generally excluded by our EBV and 
KQ85 masks. 

To investigate the pro perties of these templates further, we have computed pseudo-C^ power spectra in the Celestial 
North and South where the IDDdI template uses exclusively WHAM or SHASSA data, respectively. The results are shown in 
Fig. IA2I The noise properties of the maps are difficult to quantify because many components can contribute, including stellar 
residuals that will dominate. However, the instrumental thermal noise is approximately known. For the WHAM survey, the 
noise level varies over the sky but is typically 0.02 R, whereas for SHASSA the corresponding amplitude is ~ 0.01 R. These 
values are overplotted in the figure. In the case of the lDDOl template, no additional smoothing has been applied, but for the 
IF03| map we have smoothed the data from the stated resolution of 6' to a 1° effective FWHIM. It is interesting to note the 
asymmetry in amplitudes between the two hemispheres. For comparison, we overplot power laws (with index and amplitudes 
chosen by eye to roughly match) as smoothed with a Gaussian beam of FWHM 1°. Both templates appear to have a roughly 
1° resolution in the South, where both use SHASSA data. The .DDD. template in the North, where only WHAM data are used, 
is also roughly consistent with 1° resolution, although de viatio ns are seen, as perhaps might be expected due to WHAM^s 
tophat rather than Gaussian beam profile. However, the iFOj template (black) does not appear to match the intended 1° 
resolution in this region, where it uses a mix of both WHA M an d VTSS data. Rather, the resolution is apparently lower, 
roughly following a beam of 1.3° instead. This implies that the lFOSi template is inconsistent with the claimed 6' resolution and 
instead exhibits more complex properties over a range of angular scales. This would appear to demonstrate the difficulties in 
combining and using inhomogeneous data sets. 

Of central importance to this paper, however, is the impact of the observed differences in the templates on our template 
fits. We have, therefore, performed analyses of a set of 1000 simulations, where we use one template to generate a sky with 
CMB plus all foregrounds, and then we fit both the correct H^ template and the alternate template to the data. In many 
regions the cross-correlation coefficients are the same, but in equally many regions systematic differences arise due to even 
small diffe rences in the templates. In some cases, it is clear whe re th e difference in results likely originates; in reg ion 4 taken 
from iDOd for example, there is an obvious bright spot in the IF03| template that is much fainter in the IDDDI relative to 
the surrounding structure (see Fig. IA3|) . In other cases, the cause of the difference in results is not apparent among low-level 
template differences. Moreover, we have found such behaviour not only in regions where the data used in the two H^ templates 
are different (e.g. VTSS versus SHASSA) but also in regions where the data should be consistent with each other (e.g. the 
common sky area covered by the SHASSA survey), implying that even relatively small differences in the sky map processing 
can affect the results in such small regions. 

With regards to the data, we find that the mean tem plate fit coefficient to the WMAP K-band data for the 33 regions is 
1.01 ± 0.02 for the lDDPl template and 8.06 ± 0.05 for the lFOSi template (in appr opriate units). Indeed, the former results by 
region are systematically lower than the latter. However, assuming that the lFOSi results are more correct simply because the 
results are closer to what we expect (based on idealised properties of the WIM such as an electron temperature of ~ 8000 K) 
is unwise. In the lDavies et al. I I2OO6I paper, partly because of the unexpectedly low cross-correlation coefficients, we applied a 
small additional smoothing to the WHAM regions. This smoothing was an attempt to make up the difference between a 1° 
top-hat and a 1° Gaussian, yet it appears from Fig. IA2l that this additional smoothing was not appropriate. 

Nevertheless, it is in smoothing that we find a solution to the issues observed here. After convolving the templates and 
WMAP data to an effective 3° resolution, then we find st atis ticall y consistent results between the two templates - for example, 
coeflicients of 9.92± 0.13 and 10.28± 0.13 for the lDDDl and lFOSi templates, respectively. 
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Finkbeiner Coefficients Finkbeiner Coefficients 

Figure A4. Comparison of the free-free fit coefficients determined between the WMAP K-band data and the lDDPl and lFoj templates 
at 1° (left panel) and 3° (right panel) resolution. 
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Figure A5. Correlation between the lDDPl and lFOSi template maps at 1° and 3° resolution for region 15 (left panels) or the EBV maslt 
(right panels). The thick blue line represents the best fit regression with slope M and intercept C for the number of pixels N. Clearly the 
scatter of the correlation plot decreases with 3° smoothing. 

The difference in results between the two templates at 1° resolution seems to be particularly difficult to reconcile with 
the power spectrum plots. One would at least expect good agreement in the regions where both templates use only WHAM 
or SHASSA data. Figs IA4| fc IA5I are therefore quite informative in this respect. It should be clear that the two templates 
are generally very consistent, although there are more outliers at 1° resolution. We therefore postulate that it is the nature 
of the estimator a = d/t^ M~^t (see Section |3| that amplifies the differences in the templates, as the denominator is 
sensitive to the square of the template. Though the numerator is roughly the same using both templates, the denominator 
differs significantly and drives the difference in a. At 3° resolution, the differences are suppressed sufficiently to allow excellent 
consistency between the results derived with either template. 

Though we cannot make a case for the superiority of one template over the other, we choose to utilise the IDDdI template 
since it was explicitly constructed to have a uniform resolution over the sky, consistent with the analysis above. Moreover, 
we wish to emphasise that, in neither case do we consider that the templates are invalid as proxies for the Galactic free- free 
emission in their higher resolution forms, but that care must be taken as to the properties of the data in local regions. This 
is particularly important given the subtle interaction between the structure in the data and type estimators as employed 
in our analysis. 
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APPENDIX B: CONSTRUCTING A REGIONAL DEFINITION MASK 

The lDavies" et al.l l|2006l ) regions were selected on the basis that one of the three principal and well-known Galactic foregrounds 
(free-free, dust or synchrotron emission) was dominant in each region. The regions were selected well away from the Galactic 
plane where the foregrounds are inevitably confused. In this paper, we extend the analysis by defining regions over a larger 
fraction of the high latitude sky, and making more extensive use of the template morphology to achieve this. Our approach 
is also intended to allow a focus on localised areas of the sky where foreground spectral variations may be traced. To achieve 
this, we recogni se that the overall G alactic emission even to high latitudes can generally be well described by a cosecant 
law of emission (|Bennett et al.ll2003bl ). with the possible exception of the 408 MHz emission. Nevertheless, we wish to adopt 
a uniform processing of the templates. Therefore, by normalising the template emission by a cosecant term, we essentially 
flatten the structure and enhance local features. We impose a 5° Galactic cut on the templates given that this region will be 
excluded by our EBV and KQ85 masks in any case, and to avoid the cosecant term causing problems. We then subdivide each 
temp late on the basis of intensity, a criterion which, at least in the case of thermal dust emission (|Miville-Deschenes et al.l 
I2OO7I ). associates variations in the physical properties of the foregrounds with mean brightness. This is then the basis for 
further processing to select the individual regions. The exact algorithm to define our regions of interest is as follows. 

(I) Consider the three external foreground templates for synchrotron S (408 MHz), free-free F (Ha) and dust D (FDS8) as 
smoothed to 3° beam resolution. Denote the emission at a given pixel of component k (S,F,or D) by T^{n). 

(II) To remove the effect of the linc-of-sight depth, we apply a sine modulation to each external template, which is defined 

as, 

r.^(^)=/r^n)xsinl61 ff H>5 
^ \0 ff |6| < 5 

where '6' represents the Galactic latitude of the given position of the sky defined by n , T*" and T^.^^ represents the external 
template and the modulated template, respectively. We do not include the ±5° latitude range above and below the Galactic 
plane. Each modulated template has a characteristic behaviour which helps to define the morphology of the regions of interest. 



(HI) Sort the pixels in each template by temperature such that 

Tsin(nO < rl(n.+i) (B2) 

and then define ten intervals, each with a tenth of the ordered pixels. I.e. for each template, we define the set of pixels 

S* = {n,} s.t. O.lAf (i - 1) < i Q.liVj 

where 1 j 10 (B3) 

This effectively divides the sky into ten roughly isotemperature regions of equal total area. We can create a map, /'°(n), of 
these sets by assigning the pixels the dummy temperature value of j that labels the set, , to which the pixel belongs. These 
maps are shown in the first row of Fig. IBll for the three templates. Where the map has a value of 10, for example, the pixel 
in the sine-modulated template map has a temperature in the top 10% of the sky emission. 

(IV) This isotemperature contour map for each template has distinct features and can be used to isolate morphologically 
interesting regions. We visit each pixel and determine in which template l''{ni) is maximum and then define a new region 
map, R{n) as follows; 

R{hi) = a max{/''(7ii)} 
f 10 if max{j''} = 

where a = i 100 if max{j'''} = (B4) 
[1000 if max{j''} = 7^ 

This map is shown with various color scalings in the 2nd row of Fig. IBll Levels 10-100 correspond to synchrotron regions, 
100-1000 are dust, and 1000-10000 are H^. 

(V) The new combined map R{n) has up to ten regions for each template, and these are reduced by a factor of two into 
15 (non-contiguous) regions simply by merging regions of similar temperature. This is shown in the left-hand plot of the 
last row of Fig. IBll Lastly, the discontinuous regions with the same value label are then split, while very small regions are 
either clumped together or omitted, leaving a total of 35 new regions. The resultant map is shown in lower-right of the figure. 
Table iBll summarises some useful quantities that describe the regions on the sky. 
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Figure Bl. Definition of tfie new regions as described in i; 12.41 Top row. the sine-modulated external synchrotron, free- free and dust 
templates at 3° resolution, each rescaled into ten roughly isotemperature regions of equal area. Second row. regions of the sky where where 
each modulated template is considered morphologically dominant. Last row. the left-hand plot shows the merged map from combining 
and maps from the middle row. The central plot defines 35 regions constructed from the merged maps. These cover almost the whole sky 
excluding the ±5° close to the Galactic plane and some areas where the method resulted in many very small regions that are difficult 
to merge meaningfully. Note that two of the regions are eventually excluded from analysis when the KQ85 mask is also imposed on the 
sky. 
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Region No. of Pixels (in N^side — map) Angular Coordinates Centroid of Region 

Default Surviving EBV mask I (in degrees) b (in degrees) I, b (in degrees) 
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41°. 6 to 178° 


6 


32°. 8 to 88° 
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(95° 





62° 


2) 
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1433 


190°. 7 to 363° 
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25° .3 to 87° 
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(302° 


2 


59° 


2) 
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2508 


286° .9 to 60° 
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15° .7 to 88° 


5 
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1 


58° 


3) 
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180°. 8 to 260° 
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37° .9 to 67° 


9 


(202° 


1 


55° 


5) 
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145 


166°. 8 to 179° 


3 


35°. to 55° 


9 


(174° 


3 


45° 


2) 


6 


670 


585 


39° .4 to 9,-)° 


6 


17°. 6 to 52° 





(67° 


6 


33° 


2) 


7 


1624 


1507 


181° .4 to 251° 





6°.0 to 56° 


6 


(208° 


5 


30° 


4) 


8 


1026 


920 


239° .1 to 333° 


3 


8°. 4 to 34° 


2 


(287° 


6 


22° 


3) 


9 


775 


260 


345°. 9 to 31° 


6 


6°.0 to 41° 


8 


(7° 


7 


20° 


9) 


10 


702 


575 


128°. 7 to 179° 


3 


6°.0 to 34° 


2 


(160° 


4 


18° 


9) 


11 


696 


594 


55°. 6 to 111° 


1 


6°.0 to 35° 





(75° 


7 


18° 


2) 


12 


1093 


279 


79°. 6 to 168° 


9 


6°.0 to 33° 


5 


(121° 





14° 


2) 


13 


421 
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243°. 3 to 277° 





6°.0 to 25° 


3 
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8 
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9) 


14 
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-28°. 6 to -6° 
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5 
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8) 


15 
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8 
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(84° 


6 


-16° 


6) 


16 
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479 


9°.l to 53° 


4 


-30°. 7 to -6° 





(31° 


3 


-17° 


2) 


17 


290 


269 


308°. 7 to 358° 


6 


-24° .0 to -10° 


8 


(330° 


2 


-18° 


0) 


18 


627 


573 


109°. 7 to 149° 


1 


-30°. 7 to -9° 





(130° 





-18° 


9) 


19 


330 


275 


229°. 9 to 268° 


7 


-30°. 7 to -10° 


2 


(238° 


8 


-22° 


5) 


20 


286 


135 


181°. 4 to 192° 


7 


-39°. 6 to -12° 


6 


(186° 


2 


-26° 


1) 


21 


899 


703 


276°. 3 to 319° 


2 


-44°. 2 to -10° 


8 


(297° 


7 


-27° 


5) 


22 


392 


354 


87°. 2 to 111° 


1 


-41°. 8 to -6° 
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2 


-29° 


2) 


23 
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496 


148° .4 to 179° 


3 


-45° .8 to -12° 


6 


(166° 


8 


-28° 


6) 


24 
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812 
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5 


-51°. 3 to -12° 


6 
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-34° 


4) 


25 
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54°. 5 to 77° 


9 


-49°. 7 to -34° 


2 


(66° 





-43° 


5) 


26 
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235° .2 to 272° 


2 


-52° .0 to -34° 


2 


(255° 


2 


-44° 


4) 


27 
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348 


197°. 4 to 239° 


1 
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1 


-46° 


7) 


28 
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87°. 3 to 123° 


9 


-55°. 1 to -42° 


6 


(106° 


5 


-60° 


2) 


29 
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7°. 9 to 35° 


8 


-63°. 4 to -50° 


5 


(19° 


8 


-66° 


8) 


30 


195 


182 


168°. 1 to 179° 


1 


-64° .9 to -48° 


1 


(168° 


7 


-56° 


9) 


31 


579 
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181°. to 305° 


2 


-67°. 2 to -38° 


7 


(250° 


2 


-59° 


9) 


.'-!2 


3!)88 
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181° .2 to 148° 


fi 
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(.3.-,-° 
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4) 



Table Bl. Definition of regions for the template fit analysis. The number of pixels contained in the region by default, and after application 
of the EBV mask are shown, together with the range of latitude and longitude values spanning a given region, and the angular coordinates 
of the region centroid. 
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3 
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.0 
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±5 
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.3 
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7 
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±1.7 




±1 


.7 
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±1 
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.7 
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±0.6 




±0. 


6 


0.5 
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5 




±2.5 


5 2 


±2. 


4 


4.6 


±2. 


3 


3.9 


±2. 


1 


11 


13.9 


±1 


.4 


3.2 


±1.4 


0.8 


±1. 


4 


-0.7 


±1 


.3 


-1.2 


±1 


.2 


12 




±2. 


1 


3.2 


±2.1 


2.2 


±2. 





1.1 


±1. 


9 


0.8 


±1. 


7 




2 6 


±1. 


9 


-0.5 


±1.9 


" ■ 


±1 


.9 


-1.2 


±1 


.8 


-1.3 


±1 


.6 


14 


■ 


±2 


.1 


-3.2 


±2.1 




±2 


.1 


-2.2 


±2 


.0 


-1.6 


±1 


.8 


15 


11.0 


±1 


.8 


2.1 


±1.7 


-0.3 


±1 


.7 


-1.8 


±1 


.6 


-2.2 


±1 


.5 


16 


6.3 


±0. 


7 


2.1 


±0.7 




±0. 


7 


0.3 


±0. 


7 


0.1 


±0. 


6 


17 




±1. 


4 


1.4 


±1.4 


1 


±1. 


4 


-1.1 


±1 


.3 


-1.1 


±1 


.2 


18 


8 5 


±2. 


2 


4.0 




3 


±2. 


2 


1.7 


±2. 


1 


1.5 


± 1 . 


8 


19 


2.3 


±5. 





-0.2 


±4.9 


-0.9 


±4 


.8 


-2.7 


±4 


.7 


-2.3 


±4 


.2 


20 


4.9 


±4. 


2 


2.1 


±4.1 


2.0 


±4. 


1 


1.1 


±3. 


9 


0.7 


±3. 


5 


21 


6.2 


±1. 


9 


1.9 


±1.8 


1.0 


±1. 


8 


0.2 


±1. 


7 


-0.0 


±1 


.5 


22 


6.1 


±2. 


7 


2.9 


±2.7 


2.3 


±2. 


6 


1.6 


±2. 




1.2 


±2. 


3 


23 


2.6 


±2. 


9 


-1.2 


±2.8 


-2.1 


±2 


.8 


-2.6 


±2 


.7 


-2.7 


±2 


.4 


24 


12.8 


±2 




6.8 


±2.4 


5.0 


±2. 


4 


3.5 


±2. 


3 


2.6 


±2. 


1 


25 


2.0 


±3. 


6 


-0.3 


±3.5 


-0.5 


±3 


.5 


-1.0 


±3 


.4 


-0.9 


±3 


.0 


26 


12.0 


±5 


.7 


5.4 


±5.7 


3.9 


±5. 


4 


2.9 


±5. 


3 


2.6 


±4. 


8 


27 


14.5 


±5 


.3 


7.8 


±5.2 


6.3 


±5. 


1 


4.3 


±4. 


9 


3.8 


±4. 


4 


28 


2.8 


±6. 


3 


-1.6 


±6.2 


-2.6 


±6 


.1 


-4.0 


±5 


.9 


-3.5 


±5 


.3 


29 


13.0 


±6 


.3 


7.0 


±6.2 


5.7 


±6. 


1 


4.2 


±5. 


9 


3.1 


±5. 


3 


30 


2.3 


±5. 


9 


-1.3 


±5.8 


-2.1 


±5 


.7 


-2.4 


±5 


.4 


-2.4 


±4 


.8 


31 


9.7 


±3. 


8 


0.8 


±3.8 


-1.5 


±3 


.7 


-2.8 


±3 


.6 


-2.8 


±3 


.2 


32 


7.1 


±1. 


1 


2.4 


±1.1 


1.3 


±1. 


1 


0.3 


±1. 





0.2 


±0. 


9 


33 


8.9 


±3. 


1 


5.1 


±3.1 


4.0 


±3. 





3.1 


±2. 


9 


2.4 


±2. 


6 


Average 


6.92 


±0. 


26 


2.11 


±0.26 


1.00 


±0. 


25 


0.15 


±0. 


24 


-0.08 


±0 


.22 



Table CI. Template fit coefliicients between the WMAP data and the lHaslam et all lll982h 408 MHz data, used as a proxy for Galactic 
emission due to synchrotron radiation, in units of /iK KJpg mHz' Results are provided for the global regions defined in Section 12.31 and 
for the 33 sky regions defined in Section 12.41 Monopole and dipole terms are also fitted simultaneously. For the global fits, NFS - North 
Folar Spur, GN - Galactic North, GNreduced - Galactic North with the NFS removed, EN - Ecliptic North, ENreduced - Ecliptic North 
with the NFS removed, GS — Galactic South, ES — Ecliptic South. The average of the 33 regions is also tabulated for convenience. 

APPENDIX C: FULL FIT RESULTS 

Tables ICll IC2I & IC3I list the solutions for all 33 regions, for the template fit coefBcients of the three Galactic foreground 
components of synchrotron, free- free and dust emission respectively. Discussions of these basic results can be found in Section|4j 
and model fits to the coefficients are detailed and discussed in Section [5l 
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Free- Free 


Region K 


Ka 


Q 


V 


W 




BBV 9.80 ±010 




4.78 ±"-"^ 




3.11 ±009 




1.29 ±0-09 




0.48 ^"■"'^ 


KQ85 10.43 ±0 20 


5.14 ±"'20 


3.40 ±010 


1.45 ±0-10 


0.62 ±"1^ 


NPS 20.99 ±12.11 


g3^±11.94 


10.35 ±11-74 


7.14 ±11-33 


5.29 ±l"-l'^ 


GN 10.31 ±0-l» 


4.99 ±0-18 


3.20 ±0-17 


1.16 ±0-17 


0.27 ±0-1^ 


GN,„duccd 10.27 ±0-1** 


4.97 ±"-1^ 


3.18 ±0-1^ 


1.15 ±0-1^^ 


0.26 ±0-1^ 


EN 13.22 ±0-26 


6.52 ±0-26 


4.32 ±0-25 


1.76 ±0-24 


0.59 ±0'22 


EN„duccd 13.12 ±0-26 


6.48 ±'^■26 


4.29 ±0-25 


1.75 ±0-24 


0.58 ±"-22 


GS 9.81 ±012 


4.79 ±"-12 


3.13 ±012 


1.37 ±011 


0.57 ±"1° 


ES 9.22 ±011 


4.47 ±*^-ll 


2.89 ±010 


1.20 ±010 


0.44 ±9.09 


± -iD.4 


1 E /I ±15.4 

-io.4 


1 r 1 ±15.0 
-io. i 


1 A a ±14.7 

-14. y 


ion it 13.3 

-io.y 


Z ii.z 


7.2 ±^-3 


K 7 ±5.2 
b. t 


A 1 ±30 
4.3 


o o ±4.4 
3.3 


3 11.8 ^'''■'^ 


1 1 ±7.6 
4.1 


o Q ±7.5 
2.0 


n A ±7.2 
-0.4 


1 A ±6.5 
-1.4 


4 17.9 


8.2 ±23-5 


y .7 ±23.0 
3.7 


1 o ±22.2 
-1.2 


^ Q ±19.9 
-4.3 


5 o2.6 


e-r, Q ±44.4 
63.8 


fir . ±43.2 

o5.4 


Rr n ±42.3 
65.0 


60.1 ±3S-2 


a IT'/ ±9-4 


o n ±9.3 

y.u 


R R ±9 
6. 6 


c- n ±3-7 

O.O 


4.Z 


7 9.4 =1=2.1 


o -7 ±2.1 


n 7 ±2.1 
0. ^ 


1 A ±2.1 
-1.4 


o n ±1.9 


8 D.9 


O Q ±1.5 

3.8 


o o ±1.5 
2.2 


n R ±1.5 
0.6 


n 1 ±1.3 
-0.1 


n n n i . 3 

y y .y 


E o ±0.3 

O.O 


o 7 ±0.3 


1 ,- ±0.3 
i.o 


n /I ±0-3 
U.4 


iU -i.o 


A ±5.4 
-4.5 


R o ±3.3 
-6.2 


R ^ ±5.2 
-O.O 


e Q ±4.6 
-D.O 


11 25.0 ='='^-^ 


T 1 1 ±0.8 
11.1 


•7 n ±0.8 
7.0 


Q Q ±0.8 
3.3 


1 n ±0.7 
1.9 


ii (ivU^oOj 15. o 


-7 1 ±1.3 
f . 1 


A /I ±1-3 
4.4 


1 fi ±1-3 
i.o 


l.U 


Lz y .y 


Q Q ±2.4 

3.3 


1 n ±2.4 


^ o ±2.3 
- 1.3 


1 a ±2.1 

-i.y 


lo y.o 


.1 Q ±0.4 
4.3 


o 7 ±0.4 
2. / 


-1 -1 ±0.4 
1 . 1 


n A ±0.3 
U.4 


14 9.8 ='='^-^ 


A Q ±0.2 
4.8 


■3 o ±0.2 
3.2 


1 A ±0.2 
1.4 


n a ±0.1 
0.6 


15 12.3 ='='^-'^ 


o ±0.7 
D.Z 


. Q ±0.7 
4.3 


o 1 ±0.7 
2.1 


1 o ±0.6 
1.2 


1 c T o zh 5 . 2 
Id I .o 


Q O ±5.1 

3.0 


o Q ±5.0 
2.3 


n ^ ±4.8 
U.3 


no ±4.3 
\.).Z 


If - U . 4 


n « ±S.9 
O.O 


0.7 ±3-3 


U.4 


n q ±5.0 
-0.3 


io iooit2.0 
lo lo.o 


o 1 ±2.0 
o. 1 


r c ±2.0 

5.8 


q R ±1.9 
3.6 


o o ±1.7 
Z.Z 


19 5.7 ±2.6 


2.9 ±2-6 


1.4 ±2- = 


0.0 ±2-^ 


-0.4 ±2-2 


20 8.0 ±"-^ 


4.4 ±0.9 


3.6 ±0-0 


2.2 ±0-9 


1.5 ±0-8 


21 13.6 ±^-2 


8.9 ±^1 


7.6 ±3-0 


5.3 ±2-9 


3.8 ±2-0 


22 7.8 ±^-^ 


1.6 ±^-'' 


0.4 ±^-0 


-1.5 ±^-■1 




23 10.8 ^^-"^ 


4.7 ±1-0 


2.8 ±1-0 


0.6 ±1-0 


-0^4 ±1-1 


24 7.4 ±"-4 


3.7 ±0 '1 


2.4 ±0-4 


1.0 ±0-^1 


0.5 ±0-4 


25 2.4 ±13.2 


1 2 ±130 


I I ±12.8 


-1.6 ±12-3 


_3 J ±11.0 


26 11.7 ±22.5 


0.1 ±221 


-3^4 ±21-3 


-8.2 ±20-8 


-5^4 ±18-^ 


27 -4.8 ±^-2 


_y g ±6.1 


-8.5 ±0-0 


-9.2 ±^-3 


-9.1 ±5-2 


28 -25.7 


-30.0 ±"3 


-30.8 ±*3-0 


-30.3 ±"1-** 


-33.9 ±37-7 


29 -42.4 ±23.4 


-40.6 ±23 


-35.0 ±22-6 


-35 1 ±22.2 


-34.3 ±20-2 


30 36.6 ±^'^-^ 


27.6 ±-l**-2 


23.3 ±-*''-2 


18.0 i-^*^-^ 


10.3 ±"2-4 


31 -19.0 ±l'='-^ 


-27.9 ±1^-1 


-36.1 ±10-^ 


-33.0 ^^^-^ 


-27.3 ±"--1 


32 9.3 ^^-"^ 


5.9 ±3 


4.4 ±2-9 


3.1 ±2.9 


2.4 ±2-0 


33 3.5 ±1^-^ 


-1.2 ±"-^ 


-1.3 ±1"- = 


-3.0 ±13-^ 


-3.5 ±12-4 


Average 9.92 ±013 


4.88 ±013 


3.23 ±0-13 


1.42 ±"-12 


0.58 ±0-11 


(a) 9.54 ±°-l^ 


4.72 ±013 


3.14 ±0-13 


1.37 ±"-12 


0.55 ±0-11 


(b) 9.60 ±"-13 


4 74 ±0.13 


3.15 ±0-13 


1.38 ±"-12 


0.56 ±0-11 



Table C2. Template fit coefficients between the WMAP data and! DDD l H(j data, used as a proxy for Galactic emission due to free-free 
radiation, in units of fiK R^^. Details as per Table ICTl Note that, as described in Section 14.31 region 11 shows evidence of anomalous 
behaviour in terms of its correlation with the Hq template. We therefore additionally show the corresponding result when the KQ85 
mask is imposed instead of the EBV mask. The average of the 33 regions is also tabulated for convenience, (a) is a variant on this in 
which we completely exclude the region 11 contribution; (b) includes the region 11 (KQ85) values instead of the standard ones. 



Foreground Analysis of 7-year WMAP data 



Dust 



R.cgioii 


K 




Ka 






Q 






V 


W 


EBV 


9.84 ±" 1« 


3.80 


±0.10 


2.16 


±0 


10 


1.07 


±0.09 


1.24 


±0.08 


KQ85 


8.23 ™ 


3.15 


±0.09 


1.85 


±0 


09 


1.06 


±0.08 


1.31 


±0.08 


NFS 


8.58 ±1" 


3.16 


±1.54 


1.75 


±1 


52 


1.02 


±1.46 


1.61 


±1.30 


GN 


9.61 ±016 


3.57 


±0.16 


1.97 


±0 


15 


0.97 


±0.15 


1.20 


±0.13 


r c d. u c c d. 


9.55 ±°16 


3.55 


±0.16 


1.96 


±0 


15 


0.96 


±0.15 


1.20 


±0.13 


EN 


8.97 ±°15 


3.36 


±0.15 


1.88 


±0 


15 


0.96 


±0.14 


1.18 


±0.13 




8.85 ±015 


3.30 


±0.15 


l.f 


i6 


±0 


15 


0.95 


±0.14 


1.17 


±0.13 


GS 


10.44 ±" " 


4.14 


±0.14 


2.40 


±0 


14 


1.20 


±0.13 


1 .29 


±0.12 


ES 


10.53 ±" " 


4.12 


±0.13 


2.34 


±0 


13 


1.13 


±0.13 


1.26 


±0.11 


1 


10.1 ±3'' 


4, 


.2 


±3.4 


2 


,7 


±3 


3 


1.2 


±3.2 


1.6 


±2.9 


2 


8.8 ±1-** 


2 


.2 


±1.8 


0, 


,5 


±1 


8 


-0.4 


±1.7 


-0.2 


±1.5 


3 


9.5 ±11 


3, 


.9 


±1.0 


2 


,5 


±1 





1.9 


±1.0 


2.2 


±0.9 


4 


6.6 ±5^ 


3, 


.6 


±5.2 


3, 


.2 


±5 


1 


2.9 


±4.9 


3.2 


±4.3 


5 


5.1 iS*^ 


1 


.9 


±9.1 


0, 


.5 


±8 


9 


0.7 


±8.7 


1.4 


±7.8 


6 


8.1 6 


2 


.9 


±2.6 


1 


.4 


±2 


5 


0.7 


±2.4 


1.0 


±2.2 


7 


10.5 ±6-^ 


4, 


.3 


±0.7 


2, 


,6 


±0 


7 


1.3 


±0.7 


1.5 


±0.6 


8 


7.3 ±''-^ 




.8 


±0.7 


1 


,8 


±0 


7 


1.3 


±0.6 


1.5 


±0.6 


9 


14.0 ±°-S 


5, 


.6 


±0.5 


3, 


,2 


±0 


5 


1.5 


±0.5 


1.5 


±0.5 


10 


y ;^ ±0.8 


2 


.4 


±0.7 


1 


,2 


±0 


7 


0.5 


±0.7 


0.8 


±0.6 


11 


4^5 ±"-^ 


2 


.1 


±0.7 


1 


,2 


±0 


6 


0.7 


±0.6 


1.0 


±0.6 


12 


5.8 ±"■'1 


1 


.6 


±0.4 


0, 


,6 


±0 


4 


0.1 


±0.4 


0.5 


±0.4 


13 


6.1 ±"■8 


2 


.1 


±0.7 


1 


.1 


±0 


7 


0.5 


±0.7 


0.9 


±0.6 


14 


10.4 ±0-5 


4, 


.3 


±0.5 


2, 


,6 


±0 


5 


1.4 


±0.5 


1.5 


±0.4 


15 


8.9 ±0-6 


3, 


.9 


±0.6 


2 


,5 


±0 


6 


1.7 


±0.6 


1.9 


±0.5 


16 


7.0 ±''''' 


1 


.4 


±0.7 


0, 


,2 


±0 


7 


-0.4 


±0.7 


-0.0 


±0.6 


17 


7.2 ±2 


2 


.0 


±2.0 


0, 


,5 


±2 





-0.2 


±1.9 


0.4 


±1.7 


18 


7.6 ±10 


2 


.4 


±1.0 


1 


,0 


±1 





0.4 


±0.9 


0.7 




19 


10.1 ±21 


5 


.1 


±2.1 


3, 


,9 


±2 


1 


3.5 


±2.0 


3.3 


±1.8 


20 


10.6 ±0-8 


3, 


.9 


±0.8 


2 


,0 


±0 


8 


0.9 


±0.8 


1.1 


±0.7 




7 c ±0.5 


2 


.7 


±0.5 


1 


.5 


±0 


5 


0.8 


±0.4 


1.1 


±0.4 


22 


9.4 ±0-9 


4, 


.2 


±0.9 


2, 


,6 


±0 


9 


1.7 


±0.9 


2.0 


±0.8 


23 


9.6 ±0-l 


4, 


.0 


±0.4 


2 


,4 


±0 


3 


1.3 


±0.3 


1.4 


±0.3 


24 


8.4 ±0-6 


2 


.4 


±0.6 


0, 


,9 


±0 


6 


-0.0 


±0.6 


0.2 


±0.5 


25 


4.9 ±2 2 


0, 


.4 


±2.1 


-0.5 


±2 


.1 


-1.2 


±2.0 


-0.9 


±1.8 


26 


4.4 ±'■■1 


3, 


.8 


±7.0 


3.4 


±6 


7 


3.6 


±6.7 


2.9 


±6.0 


27 


12.5 ±-l 'l 


6, 


.3 


±4.3 


4, 


.8 


±4.2 


3.9 


±4.1 


3.8 


±3.7 


28 


9.2 ±2 


4, 


.4 


±2.0 


3.3 


±2 





2.9 


±1.9 


2.8 


±1.7 


29 


6.3 ±13 '' 


-1. 


1 


±13.2 


-5.3 


±13.0 


-5.3 


±12.6 


-5.3 


±11.3 


30 


10.5 ±3-5 


5, 


.0 


±3.5 


3, 


,9 


±3 


4 


2.8 


±3.3 


2.8 


±2.9 


31 


12.9 ±3 6 


6, 


.4 


±3.6 


4, 


,8 


±3 


5 


3.5 


±3.4 


3.4 


±3.0 


32 


8.8 ±0 


3, 


.1 


±0.9 


1.9 


±0.9 


0.9 


±0.9 


1.0 


±0.8 


33 


9.2 ±2 1 


3, 


.3 


±2.0 


1.5 


±2 





0.6 


±1.9 


1.0 


±1.7 


Average 


8.62 ±0-11 


3.24 


±0.14 


1.81 


±0 


14 


0.93 


±0.13 


1.13 


±0.12 



Table C3. Template fit coeffieients between the WMAP data and the lFDSl model 8 data, used as a proxy for Galactic emission due 
dust radiative processes, in units of fiK /^Kp^gg. Details as per Table [CT] 



